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ABSTRACT 


Maintenance of continuous electrical service by a 
power company is possible only if its power system is stable 
when subjected to any disturbances which change power flow. 
Extensive study is necessary to guarantee this continuity of 
Supply. The study is termed power system stability analysis. 
The analysis, in turn, is wholly dependent on a non-linear 
differential equation - the swing equation - which describes 
the motion of each synchronous generator on a power system. 

The advent of the digital computer has enabled power 
systems engineers to adopt new and improved methods of system 
stability analysis. Runge-Kutta integration methods have 
replaced the step by step methods previously used to solve 
the swing equation numerically. 

Two Runge-Kutta methods are in popular use in this 
application: 


1) the method of Johnson and wara (2), 


2) the method of cili!t®), 

The intent of the writer is to analyze these and various other 
integration techniques with a view to finding a procedure which 
gives optimum accuracy of solution of the swing equation, with 
minimal computer computation time. The writer has derived a 


new Runge-Kutta scheme to achieve that objective, and this 


method is discussed in this thesis. 
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I. INTRODUCTION 


1.1°-The Swing Equation and Its Solution 


The swing equation is a non-linear, second-order dif- 
ferential equation derived from the fundamental laws of mech- 
anics +) , It describes the motion of each synchronous machine 
on a power system. Parameters needed for a description of the 
basic form of the equation are as follows: 

Py is accelerating power, the difference between net mechan- 
Toe rTnpuec Ai and electrical power output Poe both cor= 
rected for losses; 

M is the"inertia constant" of a machine, a term proportional 
to one representing kinetic energy at rated speed; and 

6 is the displacement angle of the machine rotor with res- 
pect to an imaginary reference axis rotating at synchronous 
speed; and 

a is the elapsed time measured from the beginning of a 


disturbance. 


The swing equation is of the form 


Mechanical input power is usually assumed constant, because 
of the relative slowness of governor action. In the most 
general case, assuming no damping and constant internal vol- 
tage, 


P. = £(6). 
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If damping and field decrement are taken into account, the 


swing equation takes the form 


In power systems design, swing curves, loci of 6 vs. 
time governed by the swing equation, are calculated to evalu- 
ate system response following an assumed disturbance in the 
generator power outputs. An analytical method of solving the 
Swing equation has not been found to date. The most popular 
procedure for obtaining the swing curves has been to obtain 
approximate solutions of the swing equation by numerical in- 
tegration. 

The most important integration methods used have been 
the following: 

a) A step-by-step method generally used in conjunction with 
the network analyzer‘)? = 
i) the acceleration, as calculated at the beginning of 
a particular time interval, remains constant from the 
centre of the preceding interval to the centre of the 
interval in question; 
jiacesthe angular velocity, as computed at the centre of 
an interval, remains constant over that interval. 
b) The number series method, recently proposed by Rao and 
Bao”). This method is based on a closed-loop control 
system representation of the swing equation; the response 


is evaluated by trapezoidal integration of the convolution 
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integral. 

c) The Runge-Kutta fourth-order method and its variations. 
The use of Runge-Kutta in this application was introduced 
by Johnson and wara ‘2? , 
The Johnson-Ward method was adjudged superior to the 
commonly-used step-by-step method by the following claims: 
a) it was more accurate than the priorly-used methods and 
took less computational time on a digital machine to ac- 
complish this accuracy; 

b) it was self-starting - every integration step was begun 
anew without reliance on any values preceding those at 
the beginning of each iteration. 


More recently ‘?"4) | 


another computational procedure has been 
used in place of the Johnson-Ward method. This is Gill's vari- 
ation of Runge-Kutta - applied to the swing equation by using 
a substitution of variable which transformed the equation into 
two first-order equations. 

Neither of these methods gives optimum accuracy and 
minimum computational time for stability studies on a digital 
computer. The writer derives and analyzes these and other 


methods, with a view to finding which is the most suitable in 


this regard. 
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II. DERIVATION OF RUNGE-KUTTA FOURTH-ORDER METHODS 


TO SOLVE THE DIFFERENTIAL EQUATION 


— = £(x,y); y(x,) = ee 


The derivation of first-order Runge-Kutta algorithms 
is certainly not new. It is presented here because a basic 
understanding of both the mechanics of the derivation, and the 
form of the notation used, is essential for a full comprehen- 
sion of the sections following. 

The general purpose of a Runge-Kutta method is to ob- 
tain an approximate numerical solution of a differential equa- 


tion. In contrast to a formal Taylor Series solution, 


hn? 
y(x + Afy= y(x®)*s=,7rh beng + — Yixe) 
a 
h? 
Fee BY cal at 5)? ie tats (2.1) 


which requires the evaluation of numerous derivatives at the 
point Xo the Runge-Kutta algorithm uses only one derivative 
of the function, and obtains y (xX, + h) at the expense of sev- 
eral evaluations of this derivative. 

The aim of the derivation is to express y (xX, +h) - 
y (xX) as a linear combination of values of the function be- 
tween Xo and x, + h.To this end we define $x.8 y i =41,*4 such 


that 


u,k, + truncation error (232) 


emitiznopis stsiuA-spivA tebxK 
otesd 5 seusosd sxed | Be 
sit bas okdevizeb sdt 16 ac 
~asdezqmoo rip2 6 02 sisneae 


-do ot at pods situA-spoud s 


supe Isitnstettib. & to inte agian 


eiiens eei1ee rolysT “tomot 5 of 3 
$ 


fl 
wa ee = igi e 


(£.8) roe + (ome et 


- 
mom 
® 


odd 46 eevisaviaed avoxemun Xo. | 
evissvined sac yimo eeen mditbxopi 
~vee to sansqxe ent a5 (al + sd 


mao 


xi ey Feo tyr 
- (a + ae meoudie 82 8b ie . tee 


where the $u,8 are constants minimizing the truncation error. 


In a fourth-order Runge-Kutta method, the term 


must correspond to the true Taylor expansion (2.1) in all terms 


up to and including those of order Hes 


“ 5 h h° (3) 
sea pa keer Ch) teh yi(x oY Nauue oe meee 7 y (x,) 
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- —hey (x0) BF oa (273) 
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To avoid the evaluation of explicit higher derivatives 


in (2.1), we set 


k, = h-£(x,,y,) (2.4) 
ky = heft (x) + ah, Hes + Bk, ) rey) 
Kk = heft (x, + ash, vn + Bik) + YK) (2.6) 
ky, = hef(x, + ah, yy + Boky + Yoko + 89k3) (2.7) 


The problem now is to find suitable values of the constants 


fa,3, Re $743 ; $53 and £3 such that (2.3) is satisfied. 
Using the notation 
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In order to correlate (2.9) with (2.2) it is necessary to 
expand the expressions for the $x into their Taylor expan- 
sions about the point (X5rY5) + 


Using the Taylor series for two variables, namely 
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Obtaining these equations is the most tedious part of the 
derivation. All that remains is determining values of our 


arbitrary constants fa,8 , $5.3 ; $18 : $5.3 and a3 such 


that 
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kK, + ugk, 


Wyk, + Hok 


aed se 
corresponds exactly to (2.9) up to and including fourth-order 
terms inh. 


Comparison of (2.11) - (2.14) with (2.9) gives the 


following relationship between the constants: 
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(1951) and set forth in Martin 


is the one most commonly used in engineering applications at 


present, 


Lc 


for solution of first-order differential equations. 
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represents y (x, TT weeW LCN an error proportional to h?, al- 


though this error is small. This is Gill's method. As an 
example of the application of (2.23) and (2.24), the above 
equations will be analyzed to prove fourth-order accuracy: 
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0 1/2 1/2 1 (1 WHE | = | 172 
V2 
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V2 
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v2 _ _ 
0 1/8 1/8 1 1/4 
1 1 
0 0 = (Gh rey ee 1/12 
: iD 
0 0 0 1/4 1/24 
0 0 F(l-+) 1/2 1/8 
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Since (2.24) is also satisfied, this verifies that Gill's pro- 
cess is a valid fourth-order process. However, this analysis 
does nothing to indicate how accurate a fourth-order method 

it is - no information is given which enables a comparison of 
the accuracy with that of other procedures. The latter aspect 


will be treated later. 
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III. SOLUTION OF THE SWING EQUATION USING 
RUNGE-KUTTA FIRST-ORDER DIFFERENTIAL 


EQUATION TECHNIQUES 


The fastest and most commonly used method of solving 
second-order differential equations using Runge-Kutta is to 


break the equation into two simultaneous first-order equations 


and solve both step by step. As set forth by Hildebrana ‘+®? 
using the basic procedure of Kutta for oy = f(x,y), 
' — 
ky = hf (xX 5,+Y,) 
h k,' 
k,' = hf(x. + —, y. + —) 
Z fe) 2 fe) 9 
h k,' 
' = SS — 
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2 2 
! = ' 
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Yi = y (xX, +h). = %. eh —(ky' + 2k,' + 2k 3" + ka')yp 
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the method becomes as follows when applied to the equations 


dy 
— = F(x,y,u) and (3.2) 
ax 
du 
mt = G(R, Vyu) 2 Soe 
dx 
1 
Wer, to copy + —(k, + 2k. + 2k, + k,) and (3.4) 
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For an equation of the same form as the basic swing equation, 


say 


ey rs Ly iy Vow), (3.8) 


the simultaneous equations are 


— = b = gix,y,u)2 and| — = f(x,y,u). ° (3.9) 
dx ax 


In order to compare different Runge-Kutta methods applied here, 


in the same manner as in (3.2) to (3.7), it is expedient to 
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derive a conversion technique. Suppose we are given a first- 


order differential equation method such as 


' ante 
ky = hf (x ,+y,) 
ck," BE etx. Ped h Ae 
2 O° 1 Eee foe ecard 
a b 
bs 2] 24 
k! = hf (x, + oh nie + — kj! + — k,') 
2 2 
a b 
7 31 1 . ap : 
k,' = hf (x, + ah ee + Sees + —— k, + — k, ), 
2 2 2 
and 
u u u u 
a 2 3 4 ; 
Vite GYa 4 ee oe k, ' tied k,! + ray hy 
2 2 
where 
dy 
y. = y(xX jeu y. =@.y.(x_ + h),,, and — = f(x,y). (3.10) 
re) e) 1 fe) Ae 


Expressed in tabular form, this is: 
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Using the generalized method of 


to apply the same procedure to 


(3.11), it is now possible 


(3.9) as was done in 
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In tabular form, the procedure is as follows: 
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Given any Runge-Kutta procedure to solve a first-order dif- 


ferential equation 


— = f£(x,y); y(x_) = fe 


by comparing coefficients with those of the general procedure 
SereOmel gists. 1)) iti is possible. to use~(3.19) to obtain 


an analogous method for use in solving 


d'y 
aw] = f(x v,y ')eewath y (X,) ies and yl x) = fe 
Note that, in contrast to Hildebrand's method (3.4) - (3.7), 


only four k's need be evaluated. This amounts to a consider- 


able saving in computation time. 


3.1 Conversion of Four Basic Runge-Kutta Processes 


The most commonly used first-order differential equa- 
tion methods are outlined by Martin. These are: 

a) Kutta's Simpson's rule; 

b) Gill's process; 

c). Strachey's process; 


d) Boulton's process. 


Kutta's Simpson's rule has been adapted for second-order dif- 


ferential equations by Hildebrana‘?® , For the equation 
a“y 
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and (2.26). By comparison with Table (3.11), 
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These constants and (3.11) completely characterize the pro- 


cess. 
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Substitution of (3.22) into (3.11) gives the required result: 
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then 
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In a similar manner as before, substitution of the describing 


constants into (3.11) gives the second-order procedure: 
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Boulton's process 


For a first order differential equation Boulton's process is 


characterized by 
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The corresponding second order procedure is: 
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2.2.4Soelution ofia Sample Problem 


A comparison of the suitability of these methods to 
the swing equation will be done using the most common example, 


(1) 


one from Kimbark 


Example Sark 


A 25 MVA, 60 c/s waterwheel generator delivers 20 MW 
over a double-circuit transmission line to a'ilarge metropoli- 
tan system, which may be regarded as an infinite bus (Fig. 3.1). 
The inertia constant H is 2.76 MJ/MVA and the direct-axis tran- 
sient reactance of the generator is 0.30 p.u. The transmis- 
sion circuits have negligible resistance, and each has a react- 
ance of 0.2 p.u. on a 25 MVA base. The voltage behind the 
transient reactance of the generator is 1.03 p.u. and the vol- 
tage of the metropolitan system is 1.00 p.u. A 3-phase short 
circuit occurs at the. middle of one transmission circuit. Ob- 
tain swing curves for the generator assuming a clearing time 
of 0.60 second. 

The amplitudes of the power-angle curves obtained 


after network reduction are: 


Condition Amplitude 
Pre-fault Lee obs Us 
During fault 0.5936), 0.U. 


Post-fault 2.06 p.u. 
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Two equations result from this data, after a change in vari- 
able: 
a) valid during fault: 


a 


tr a OPE Seta Fi SOs Sa terc (3525) 
aT 


b) valid after clearing: 


2 


— = YC 30) =— 2. e0Le eine on ‘o.2) 
aT 


In each case, 
6 is the displacement angle of the machine rotor with respect 
to an imaginary reference axis rotating at synchronous 


speed; 


FIG. 3.1 


SINGLE- LINE DIAGRAM OF A POWER 
SYSTEM SHOWING REACTANCES IN 
PER- UNIT. 
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is proportional to elapsed time measured from the beginning 


of the disturbance. 


From the pre-fault .data;-it-is*found-that at T=0, 


for) 
Il 


0.316 electrical radians 


dé 


dt 
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I 
= 
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el. radians/second Coes) 


Since a damping term is not included, both equations are of 


the general form 


a*6 

— = £(6); with Solu =o 6°(0) = w 

ar? fe) fe) 
Here, £(6) = 0.855 - C sin 6, where C has one value during 


the fault and another after the fault is cleared. 


The problem is quickly and accurately solved ona 


digital computer, using the flow chart shown in Fig. 3.2. The 


four numerical integration methods compared are: 


a) 
b) 
CG} 


d) 


the Kutta procedure, (3.21), 
the Gill procedure, (3 823)y 
the Strachey procedure, (3.25), and 


the Boulton procedure, (3.27). 


Seven time increments are used in each case, varying from 


h 


= 0.20 sec. to h = .003125 sec., to give an.idea of the 


relative stability of each method. For each value of h inte- 


gration is carried out from-t'= 0~to’ t = 1.0 second. The 
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.( INPUT 
NETWORK 
PARAMETERS 


2.[SET THE INCREMENT 

h TO INITIAL VALUE 
SET TIME =0 | 
SET INITIAL CONDITIONS 
8= 8), w= wo. 


4.};PERFORM ONE 
COMPLETE RUNGE - 
KUTTA ITERATION ON 
FUNCTION 


3. 


S.} TIME = TIME +h | 


HAS CLEARING 
TIME BEEN YES 
REACHED ? 


NO 


ALTER FUNCTION 


HAS ONE COMPLETE 
CURVE BEEN 
CALCULATED ? 


YES 


IS STUDY COMPLETE? 


pd i 


8. OUTPUT 
RESULTS 


FIG. 3.2 ACCURACY COMPARISON FLOW CHART 
OF EXAMPLE 3.1 
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actual computer results are shown in Appendix (C). Figure 
(3.3) compares the curves obtained by each method for a step 
size of 0.20 seconds with the true curve. 

From the figure, it can be seen that the process of 
Strachey, (3.25), gives results which follow the true curve 
most closely. In addition, this method tends to give a more 
"pessimistic" view of the swing curve; the other methods sug- 
gest that the equation solved is more stable than, in actuality, 
it is. The pessimistic view is safer. Inspection of the di- 
rect computer results (Appendix (C)) indicates that, as h is 
decreased, the Strachey process tends to converge to the true 
value of the solution nore saickit than do the others. For 
these reasons, the process of Strachey is superior to the 


others in this application. 
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METHOD 


STRACHEY REPEATED:! 79,133 


O CORRECT CURVE: 54.283 


“TS 


RELATIVE ROTOR ANGLE OF GENERATOR, 5, elec. degrees 
> co 
© [@) 


0 0.4 | 
TIME AFTER FAULT, seconds \ 

2 BOULTON REPEATED: |-19,877 

© KUTTA REPEATED: | -26.636 

CD BILL REPEATED: - 34,415 


FIG. 3.3 SWING CURVES OF EXAMPLE 3.1 - CLEARING TIME 
0.60 SEC 
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3.3 Application of Merson's Variation to 
Second-Order Differential Equations 


(12) 


An alternative procedure suggested by Merson allows 
an estimate of the truncation error to be computed at each step. 


This is, to evaluate five k's as follows: 
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for a first-order differential equation, can lead to an over- 


alivgain of efficiency ‘+) , 


despite the additional evaluation 
of f required at each step. 


In order to adapt this process so that it is suitable 
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As a means of comparing the accuracy of the "Merson 
Repeated" procedure with methods described previously, repeat 
Example 3.1, assuming a fault clearing time of 0.6 seconds, 
using the iteration method of (3.34). 


Again, the equation to be solved is of the form 
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The flow chart used to program the problem for a digital com- 

puter is the same as the one of Fig. 3.2. Step 4 in that 


figure encompasses the above operations (3.36). 
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The Merson iteration procedure is compared with the 
Strachey procedure in Fig. 3.4 for time increments h of 0.10 
and 0.20 seconds. Using the larger value of h, the Strachey 
method is more accurate. The Merson method, however, con- 
verges more rapidly to the true solution when h is decreased. 
It must be borne in mind, when comparing the procedures, that 
the latter is less efficient - it requires five evaluations 
of the function per iteration, one more than Strachey's. This 
proves costly in the case of the swing equation. The computer 
requires many addition, multiplication, and exponentiation 
steps to calculate a sine function each time the subroutine 
is entered. 

The real advantage of the Merson method lies in the 
available estimates of the truncation errors committed at 
each step. Appendix (C) gives the error values at each step 
for the example, and some are plotted in Fig. 3.5 to give an 
idea of their magnitude and variation. Table 3.1 compares 
the results of a computer program using the interval changing 
criterion of (3.32a), as set forth by Merson (with the second 
order integration set forth in (3.34) and (3.36)), and the re- 
sults of the basic program used in Example 3.1, for Strachey's 
method. The number of entries into the function, given in 
Table 3.1, provide an estimate of the relative amount of com- 


puter time required for each. 
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IV. DERIVATION OF FOURTH-ORDER RUNGE-KUTTA METHODS 


FOR SOLUTION OF THE DIFFERENTIAL EQUATION 
dy 
75 a EAX¢¥ ey ') 3 y (x,) = y 


Peyetxs) 90 (4.0) 
ax f 0 


fe) fe) 
The classic derivation of the Runge-Kutta algorithm 
for a first-order equation was outlined briefly in Chapter 
II to indicate the manner in which the technique was originally 
formulated. The same pattern will be followed, in more detail, 
in building an analogous constraint matrix for methods di- 
rectly applicable to solution of the swing equation. 
In developing Taylor series expansions for functions 
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The Taylor series for one variable gives 
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Although the object is to derive a fourth-order method, terms 
in he will be retained to give an estimate of a good part of 
the truncation error. 

Substitution of (4.3) - (4.8) into (4.9) gives the 
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The expression for k, is carried along with this derivation 
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for purposes of finding a method similar to Merson's for 
second-order equations. 
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Multiplying through this expression by h, and recalling that 
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The $7553 and 61,33 are constants appropriately chosen such 


that (4.18) and (4.19) correspond to (4.10) and 


h>. 


(4.17) 


to order 


The $x,3 must now be expanded using the Taylor series 


for three variables set forth in (4.1) and (4.2) 
If it is assumed that 
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Collecting terms gives the final expression for ko, 
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These results may now be compiled to give the constraints for 
the second-order equation Runge-Kutta algorithm. Until needed, 


terms in k. will be deleted from the compilation. 


2 


pty fe? a? + Tg) -s eur + ate, wee 
7 ett 
"Gq" 2 oY wuiedis 
S Sear | 
= Bi] ae os ry : ae a : a - 
(ésu8). 5 Co ae RO gape + ¢38 
: z aes : : 
asda. : ¢® 
pk + go + 8) + 2°@ ea) — + Boge): nn 
45 iia Tis ad | mary) 
€ 
fa — + {iG oe »& 230)S° Ms fiend 
Har eo) — + (gbgat ete +e? 2) om 


+ -7';o + es ale # ee 23 ay 


elt gre + th + cathe + 98 


is ain. od. * Jae 
bad pete lero. t g 39) + ne ee 
2 geo 


To + 630) Paeke * yitites + 3 * 


sot ednisisenol ert evin 0+ be gm 
alt — 10% = 
baboon Lisnu ate  Btta 


« 49 « 


For correspondence between terms of (4.18) and those of (4.10) 
to an order of hn? in y (x, + h), we have: 
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Likewise, for terms in (4.19) to correspond to those in (4.17) 
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Provided the relations set forth in the upper submatrix are 


satisfied, and provided 
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the following procedure is a Runge-Kutta method of fourth- 
order accuracy in h, suitable for incremental solution of the 


Swing equation as set forth in (4.0): 
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V. SELECTION OF THE RUNGE-KUTTA SECOND-ORDER 
METHOD MOST SUITABLE FOR DIGITAL COMPUTER 


SOLUTION OF THE SWING EQUATION 


With (4.51) - (4.53) in hand, it is possible to syn- 
thesize an endless number of numerical integration techniques, 
all of which are appropriate to the general equation 


ay dy 
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ax ax conditions. 


As an illustration of the procedure to be followed, the 
technique introduced by Johnson and Ward in "The Solution of 
Power System Stability Problems by Means of Digital Com- 


(2) 


puters", will be developed. 


wee Development of the Johnson-Ward Formulae 


Maintaining the notation of Chapter IV, set 
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This is in accordance with the original procedure of Kutta: 
the function is evaluated at the beginning of the step in- 
terval, twice at the mid-point of the interval, and once at 


the end. 
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(4.51), becomes that of (5.3). 
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From this information, noting that the equations expressed 
by the upper sub-matrix must hold true to give accuracy of 
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Since the right-hand side of the upper sub-matrix in (5.11) 
corresponds to that in (4.51) the truncation in y and v is 
of order fae To obtain and estimate of the truncation error, 
the lower sub-matrices of the right-hand sides of these equa- 
tions must be examined. Let TR represent the sub-matrix from 


(4.51) and JW represent that from (5.11). Then 
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To evaluate the error of order heh first of all recall that 
if each row in JW is identical to the corresponding row in 
TR there is no truncation error of 0(h>) in the process JW 
represents. Also recall that the first column of these sub- 


matrices refers to truncation of Ysi the second column to 


truncation of a= hy,'. If we set 
stale 

a = error component due to term ID Die 
h>? 

b = error component due to term Pies See 
h> 

c = error component due to term 7 DfDf.. 
h> 

d = error component due to term 7- £ bE (5.14) 
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e = error component due to term 7-— 1 © a 
h> 
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g = error component due to term 75— Ci iens 


recalling the terms from which the matrix was derived, we can 


express the truncation error components as in (S15) and 
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From (5.16) it is apparent that the truncation error in f 
is not of order ae but of order He This is poor for our 
purposes, as shall be shown. The process requires a very 
small value of h to give satisfactory results. 

Another, more accurate procedure for second-order 
differential equations is given by Collatz ‘®) , To the author's 
knowledge this has not been previously applied to the swing 
equation, yet it is more accurate than the method set forth 
by Johnson and Ward, and requires less computing time when 
applied to a function which does not explicitly contain OC), 
such as the swing equation in its most commonly-used form 
(with no damping term). The Collatz procedure is set forth 


below: 
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The error expressions for the Collatz procedure are 
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Note that here e(y,' = yeu! = O(h?), whereas in the Johnson- 
Ward procedure e(yy' - vA), = o(h*). 
Example 5.1 


The two methods are compared by a numerical integra- 
tion of the swing equation of Example 3.1, assuming a fault 
clearing time of 0.60 second. Computer results for various 
values of step length, h, are given in Appendix (C). Figure 
5.1 shows the swing curves generated by each method for 
h = 0.20 and h = 0.10 second. 

The Collatz method is clearly the superior in this 
application; it requires just three function entries per 
step because ky = ky in this case. The other procedure, in 
addition to requiring much more computer time to solve the 


problem, is far less accurate. 
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5.2 Formulation of New Second Order Runge-Kutta 
Procedures Pattern After Common 


First Order Procedures 
A. A Method Based on Strachey, Eqn. (3.24) 
The first-order Strachey process for solution of 
— = f(x,y) 


is characterized by the following integration scheme: 
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When forming a second-order method, for solution of 
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From (5.21) we may set forth the equations that must hold for 
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where y' does not appear explicitly, the integration may then 


be accomplished in the following manner: 
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Note that only three function entries are required per itera- 
tion, as compared to the usual four. 

Inspection of the lower sub-matrices of (5.19) in- 
dicates that setting 8 = B, + 54 will minimize error com- 


ponents of order h?. 
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and the completed second-order procedure is: 
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In a manner identical to that of the previous section, we 
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this process is more accurate. 
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B. A Method Based on Gill, Eqn. (3.22) 


In a manner similar to that of the previous section, 


we increment dy/dx in accordance with the Gill procedure for 


first-order equations. 
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Error components in this procedure are as follows: 
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C. A Method Based on Boulton, Eqn. (3.26) 


The constraint matrix equation is formulated, as be- 
fore, using as a starting point the characteristics of the 


first-order Boulton scheme: 
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which necessitates, as in previous sections, the use of (4.51) 
to derive a suitable method of integration. 


The constraint equations applied to (5.35) are: 
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The three second-order methods are compared by a 
numerical integration of the swing equation of Example 3.1, 
assuming a fault clearing time of 0.60 second. Computer re- 
sults for values of step length h from 0.20 sec. to 0.003125 
sec. are given in Appendix (C). Figure 5.2 shows some of 
the resulting curves. 

As predicted, the process of (5.25) gives poorer 


results in this application than do the other two methods. 
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The method based on Gill, (5.32) reduces to the Collatz 
procedure of (5.17) since the swing equation in this example 
has no damping term. Equation (5.46) gives superior results - 
better than even the Collatz method. It is worth noting that 
each of these methods takes less computer time than did any 
of those of Chapter III. Each requires only three function 
entries per iteration, in comparison to the previous four. 
Figure (5.3) compares the three most accurate methods 
from the previous section and from Chapter III, showing the 
process of (5.46) to be the most suitable of the three in 


this application. 
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VI. INTERVAL CHANGING PROCEDURES FOR 


SECOND ORDER DIFFERENTIAL EQUATIONS 


It was observed in Chapter V that Runge-Kutta pro- 
cedures derived especially for use with second-order dif- 
ferential equations gave results far superior to those de- 
rived by reducing the swing equation to two first-order equa- 
tions and applying a first-order Runge-Kutta procedure twice. 
This observation gives an optimistic picture of the possibil- 
ity of solving the swing equation using a second-order pro- 
cedure in the same manner Merson's first-order procedure was 
used. 


To this end, an expression for k, was carried through 
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in the second-order derivation of the Runge-Kutta algorithm 


(Chapter IV). If k,. is included in the formulation of the 
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The operational table, defining the symbols, is as follows: 
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Two interval-changing procedures will be developed. The 
first has, as a basis, the first-order Merson integration 


scheme; the second is based on the Boulton procedure. 


A. A Process Based on Merson, Egn. ek fee 
The first-order Merson procedure is given by the 
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The above form was easily obtained from (3.31) by assuming 


the first-order equation 
——s = f(xXyyG)» 


If we multiply both equations of (6.4) by > and let 


kK; = 3 ki! , the {3 are evaluated in the manner 
2 
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and the integration is accomplished in the manner 


5 
hy,' = hy,' + F(k, + 4k, + kg) + O(h?); and 


15¢ (hy, ' >: BY Am ~ 2k, - 9k + 8k) - Ke 
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The equation we must solve is 


d*y 
ane = SAX YO) ' (6.8) 
ax 
so that the expressions of (6.5) must be altered somewhat to 
allow for the express appearance of y in the function. 


From (6.5), referring to table (6.2), we have 
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Substitution of these values into the constraint matrix 


equation (6.1) gives equations (6.10). 
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From these equations, the following relationships 


are determined - 
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Again, the swing equation with no damping term will be con- 


Sidered. Then 


and the terms involving fu may be neglected. In order to 
save one function entry per integration step, it is expedient 


to place the constraint 


on the system of equations. Let us define ua such that 


the iteration error in y, ely, - me? is given by 


e(y;-y.) = ; w4k, + 0(h*) 
enn 4=1 ah Fe 


6 
= e,a oo Eb 4 E34 + OCH) » (6442) 


Error components a, b, and d are identical to those used in 


the preceding two chapters. 


Our system of equations then is in 
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In the case under consideration ky = Kk so that Uo + 5 = -9 


is the only constraint on Uo and U3 The coefficient B 5 was 
chosen to be zero, giving 
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The operational table, deleting Kk and setting k3! = k 


Ky! = Key ki! = ky and k,' = kK. = K becomes with Vv; 
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B. A Process Based on Eqn. (5.46) 


(6°27) 


Using a somewhat similar technique to that of Chap- 


ter VI (A), the following operational table was derived: 
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In each case kK. Nt E(x, ,yi+y;')- (6.18) 


It should be noted that the integration procedure is identical 
with that set forth in (5.46), which proved the most success- 
ful of the methods tested. The method is more versatile than 
(6.17) because it is applicable to a function containing the 
first derivative explicitly. In the latter case, however, 


five function entries are required per iteration, since 
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Example 6.2 


Compare the interval-changing procedures of Eqns. 
(6.17) and (6.18) with respect to their speed and accuracy 
by integration of the swing equation of example 3.1, assum- 
ing a fault clearing time of 0.60 seconds. 

The criterion used to evaluate the necessity of an 
interval change during the integration was identical to that 
setimorth by Merson, as in Fqn; (3.32a). “EBadgure 6.1 shows 
the procedure derived from (5.46) to be more suitable in this 
application. It requires six less function evaluations than 
does the other, yet achieves a more correct result at time 


one second after fault. 
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FIG. 6.1 SWING CURVES OF EXAMPLE 6.1 — CLEAR. TIME 
0.60 SEC 
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VII. SUMMARY 


With the increasing use of the digital computer in 
power systems planning, engineers have for the most part aban- 
doned the step-by-step integration method formerly used in 
conjunction with the network analyzer to evaluate power sys- 
tem stability. The method of Runge-Kutta has come to the 
fore. 

Presently the most popular variation of Runge-Kutta 
is the method of Gill. The swing equation is separated into 
two simultaneous first-order differential equations, as in 
Chapter II, and Gill's method is applied to each. When using 
this method, four successive approximations are made for 6; 
each approximation dependent on the preceding one. In the 
case of a single machine problem, which arises commonly, each 
approximation requires the calculation of a sine or cosine 
function. In a multimachine case the number of sine/cosine 
entries made per integration step rises rapidly. The evalua- 
tion of these functions is relatively time consuming on a 
digital computer, because of the nature of the procedure 
which must be used - hence, it is important to minimize the 
required number of function entries, while retaining desired 
accuracy of solution. 

In Example (7.1) the "Gill Repeated" procedure is 
compared with four others, to provide an illustration of the 
rapidity of each. The method of De Vogelaere, introduced in 


this example, is set forth in Appendix A. 
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Example 1 BLE 


Using the equation of Example 3.1 with fault clearing 
time 0.60 second, perform as many integration steps as neces- 
Sary to achieve an accuracy of 1.0 el. degrees in 6 at time 
1.0 second after fault. Compare the computational time for 
each of the following: 

i) Gill Repeated procedure - Eqn. (3.23) 

ii) Johnson Ward procedure - Eqn. (5.12) 
iii) Collatz procedure - Eqn. (5.17) 
iv) Second-order procedure - Eqn. (5.46 

v) De Vogelaere procedure - Eqn. (A.1-.5) 

The problem was approached by first calculating a 
swing curve with an interval size h = 0.20 second. Another 
curve was then calculated, using the same data but with step 
size 0.10 second. The values of delta at time one second 
after fault (5, o) on each curve were then compared - if 
these values differed by 1.0 el. degree or less, the calcu- 
lation was stopped and the curves were printed out with a 
statement indicating the procedure had converged to the 
prescribed accuracy. If the values of §1 06 differed by more 
than this amount, the step size was again halved, a new curve 
was calculated, and the test again applied. 

The program was written in Fortran IV language and 
executed on an IBM 7040 digital computer. For some of the 


methods it was found that total calculation time was less 
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than one second. The entire problem was thus solved 100 
times for each method, and the calculation time divided by 
100. To give a more accurate comparison between each, Table 
7.1 sets forth the results of this study, and computer output 
is shown in Appendix C. 


Table 7.1 


Results of Example i ipent 


Method of Equation Sa Doe eee AGA ero 
Acceptable value of 
$4 0° el. degrees Tae a Oe a0 7  O UL Ly s 


NOwe function en— 
tries necessary 84 164 6 20 85 
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value encountered OS Oooo OU, 207 TAS iy Bes 0.255 


Necessary step 
size U5 UD eee 0.05 U2 0.025 


7040 execution 
time in seconds 0.68 eee OR er Paes Onis 


The method proposed by the writer takes less than fourteen 
percent of the time required by the method of Johnson and 
Ward and about twenty-five percent of the time required by 
the "Gill Repeated procedure" to solve the problem to the 
accuracy prescribed. Computational time increases rapidly 
as the number of machines involved in a stability study in- 


creases, so the time saving becomes considerable indeed for 


multimachine systems. 
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Example died 


For the system shown below, the swing equation is 


2 


1 av é 
TPE (kaye 0.20 - 1.20 sin 6 per unit. 
Pee 240i idt 


With the fault applied, this becomes 


1 ds 


TEE ee tar 0.200 per unit 

te, 2 eOredt 
In this case, a three-phase bus fault is applied on bus A at 
t = 0 and removed at t = 0.30 second. Using each of the 
methods considered in Example 7.1, calculate and plot swing 
curves up to t = 1.00 second using time increments of 0.10, 
0.05 and 0.025 seconds. 

The computer output was written directly onto tape 
and an automatic plotting subroutine was used to give the 
curves shown in Fig. 7.1 - Fig. 7.5. 

From the figures, it may be seen that the second- 
order method of (5.46) converges much more rapidly than do 


the others. The larger time increment again lessens the 


total computational time. 


E,=1.2 9.6°p.u. 


3-phase 
H,=0.883 p.u. ee 


eee ee 


$s A aud m0 peitage et | aici Solel a 
ors to dose pated eters 69" 


piiwe ‘4tolq bas etsivolso whofe . 
.0f.0 to atnemexont \omts ae wooee 00.4 = 


a bore 
Hi 


sqst odno “yidoott nossin one signe 


sdt svip oF 5seu 25W ~~ 


ob asit vibiass etom don wp 10D (at 
Avs 
efit anseeol rises snonenoni “2 


a 


= 
Oo 

<9; 
“ts 

n- 
Qa 
w 
ao 
wn 


rs ec? i Rae epee aa 
Cf TESS OF SE OSH GE #% : 


te.  Peseres Feeds sys « 
etn ees eRe , 


ee ed a Ed 


- ~ 
pf. se 
s oS Peet «=~ ere einem re 
. 4 i ELS IS Eee OE " 
* os West Se ak 
Shes. wae al ~ © 
P ca - »%: aoe oe 


4 a) oe Papecmer ys «= 


? 
‘ 


pe WF 


So 


7 sony . a 
=) ae ee 


noe 


ee 


tee 


se a 


—_- ee he 


4 
ea ~ i é tae _ ° 
a ee  } a os 
ee - ae 4 S 


33 


i 


| 


| 


- 7 . a -— ‘ * = ae ee 
fe cSe C2 BSR als POSTE tee -? 
ow Vwreerwr Cu POG st wera 


* oe on 


oe ea 


P 2) 1 emilee 
~ae: 
i St 3 


a: 


ee 


» 


i 
tf 


- 
iE a 


‘ 


ey a etd 


4 
‘ 
i 
{ 


3 ee 
Dah 
e ) he ¢ 


te até 


N 


Bnd 3 


os 


SHIMLA F 
moot 
| 
: 


neh 


ty, So 


me 


eM 


as So 
t 
| 
| 
9 
4 


' 
' Liat 1 L 
ie aa 
1 ct ean 
Nai hed 
’ : 
TAaAG 
| 


94 


tenneeenee 


| 


on ae WeEPT-S BaSee es 
> Zt Soe ees <+ oC ee ee 
SRE RSLS? Val we 
ee ee eee ee ee 
= ee Ce OTe FA 


ae te. eT 


om 


4 
S 


fe ere 
i ve 


esto 98a" Shta | 


em 


1 i 
y¢ ' 
’ ‘ 


a 


SE eg ao tg ae 
> : +0- ether 
Se 


TIPE STILL ttt tec tt tcc He -+- ++ SPSS KSTC SRP SR SF SSORSSESS Sess et 
Be : SSURSGRSCcnecnescarsnesneSGGunse se 
SOOREOREeaennoEoneaanance 
PEELE LITT ELI I HEY A HE Be 44-44 He i ams 
HEHE Ht ti 7 Se Eee 
Seugnnauan agcee HE HH He 
Aseoue ae 


aaa 
eee 
. Ll aaa LL 
ee 
7. 
i. | oo 
a ee. 
i _ a 


oo CC 


oR LLSMASE ORO UBEEEN SES 
Savee at esae Bec tennese 
5 

4 


ecu 
avon 


“TL 1.t4.5.i000 


& 
z 
Sepeveeses 
Pasa geaned g 
) a. he Oe Rae a em SARS OBES UBD 
ie ie s 


Suaeauen 
syn ee ave 
ap Pe Ae ae 
sPEC HG HAA 


va ek 
3 


» 
a 
y 

- 

s 
s 
w 
a 
> 
* 

> 
a 
= 


aa cas 


3 Gseumarues 

: % BW OLRER A A 

Su seewe a a s e wa 5 2RECSP SCBSeL THT Oe a> 

<CRAGL SE CE! ORRSCH Sou4A eee » a 2 HOVER BSA PEs Pe TAS ws 
{aeCek ores eae > 


ee 
pee 


i cee Sere 
” S cqmweae: , ree 
T dS SCUSS ATI SOS 
= 


Sparen ws 0) oe ames 


See Rd sen ead -- 


i mmw ae 
t 


wre wm Rtas - - _— + . . 2 4 
EEG O25) Eeeee eeeee: ees Fe see So as ee 

Be SE RT OH - 
taeresD 


MS Severs’ 

pes CU Pi. 

TPRS2 WSs PY SS ES 
SEPT 6a 


“ 4 oe ~ - “ - - . 
+> ars ~ ie a ae a ge > 
= : ~—s : : a 
: I cides mes re ~4 
‘ 7 ba rH ' 
< . ae oe OFS ae SS cay 
* = ees | 4 bund aeons Z a tae 
2 , ore - — . ree ‘ 
7 4 : 


TaaeSkkt olaas 


ERE 


ad see 


pe mY aoe : —< S Army 4 
— : eer gt +; nabareteaiind, Lame deacon 
Se eet a pS Lees Saee ——4 ~ stom = 
a oe ee ee el eo —s . - eS Oe 


: oS SSS 
Se a eae at —+ oe ee ow et - a 
——— een Fear TE ea ONE a —+ 4 re | ~ 
ak ee pe hee fee ee ee ee ee 
tac Std rete ce ROE 23 Ome ee TRITK 6 ERATA I MEN EEE os 
ote a iad a ee =F ee — 
er ’ -4 Pe ee . . 
a OF Pane ee ae F = - Eo se ae ee o~e ee 
‘“ > ee ' + 
, Sant: rs wate 4 bw te 
__ rd — 
= SeSSeS Pe Shae toe eles 


és 


ee See 


icsntin stray Seen te 
+ 9 a ee = se and» a —— ee a ee ee ee — oe 
1 ae : re ot et s 


im dete Shah “ae eae at ana fost a ae eh ——— ie tk “ 
: 
_ aa eS >. 


og ie ge pes ee ee ete ae 


wens moos ewmeee eee Nats nS) oso 


[ee ee et ee ee a st ee ee 
‘ 


a es 


Sea SGN oe <= oe ies ace ol 7 F292 


ea passe nals 


PRR S Ce rhsy LHS 


- 96 - 


Many Runge-Kutta procedures have been presented. Most 
give both increased efficiency and accuracy of integration 
over methods currently employed. Notable among the improved 
methods are Merson's variation, De Vogelaere's variation 
(Appendix A), the Collatz method, and the second-order tech- 
nique of (5.46). In the examples studied, the method of (5.46) 
is superior. Restated in a step-by-step form to be applied 


to the equation 
2 
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2 k 
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The nature of Runge-Kutta techniques precludes any 
guarantee that the seventy-five percent to eighty-six percent 
time saving observed using this method could be accomplished 
in all cases. The scheme, however, is likely to give consider- 


able benefit in most applications similar to those presented. 
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APPENDIX A 
APPLICATION OF DE VOGELAERE'S METHOD TO THE 


SOLUTION OF POWER SYSTEM STABILITY PROBLEMS 


(14) 


A numerical integration method devised by R. De 


Vogelaere in 1955 has achieved much popularity for solution 


of second-order differential equations such as 


az 


5 += S44; 9%, 
ax 


Ke 


in which the first derivative does not appear explicitly. 


Three relationships are required to describe the basic process- 


od 1 : 4 eee - 4 
Yn =. Yo + 5 hy,’ + Banht At tsondwone (22. (A. 1) 
= i hy! c= Be (£ gh PR) O(b>) , and (A. 2) 
Yi ~*~ Yo Yo 6 O L ae : 
ee eee dean (cP aA tte othe) (A.3) 
Yi Yo 6 O +5 1 5 
For the initial integration step, note that 
(Yo-dh, Y (%e- th)? 
fy = £(x_y,Y_y) = E(x oh, y (X,o¢h)) (A. 4) 


is required. A formula suggested by De Vogelaere which is 


sufficiently accurate for an estimate of this quantity is 


P 1 DE 3 
Sev RSE hy," t+ ght, +Olh (A.5) 


The procedure is well-suited for digital computer application. 
Only two function evaluations are needed per integration step, 
and the step size may be changed with no extra steps required. 


A flow graph illustrating the technique is shown in Fig. A.l. 
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Note that in the calculation of Ko, the value of k, from 
the previous interval is used. This is a drawback of the 
De Vogelaere procedure that Runge-Kutta techniques do not 
have - the fact that the iteration is not self-starting. 
When comparing this method with others with respect 
to accuracy and time consumption, the fact that only two 
function entries are required here, per iteration, must be 
considered. As a standard of comparison, the total number 


of function entries required for integration of an equation 


over a given range is used. Table A.1l compares De Vogelaere's 


INITIAL CONDITIONS 
veh2/2 


k, zuf (85) 
ko = uf (Sp hwg/2 +k, /A) 


k= uF(Bq thug /2 +/4ky-k, | /2) 


8 = 85 +hwo +(k, +2k, )/3 
hwy = hwy +(k,+4ko* ka )/3 eth 


CALCULATIONS COMPLETE ? 


FIG. A.) DeVogelaere's procedure to solve 


d?8/dt? = £(8) 
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method, the presently-used Gill repeated procedure, and the 
second-order method developed by the writer, Eqn. (5.46). 
The problem solved is that of example 3.1, with a fault clear- 
ing time of 0.60 second. 

The Gill Repeated procedure clearly requires more 
time than the others to do a less accurate job. The De Vogel- 
aere procedure is a definite improvement over Gill's, but is 
not equal to the second-order method of (5.46) in this applica- 
tion. (5.46) is self-starting, which De Vogelaere's is not, 
and (5.46) has provision for solution of the swing equation 
with a damping term included; De Vogelaere's has not. The 
method is a definite improvement over those used to date, how- 
ever, and no discussion of numerical solution of the swing 
equation would be complete without reference to this recently 


devised technique. 
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APPENDIX B 
METHOD OF APPLICATION OF RUNGE-KUTTA METHODS 


TO A MULTIMACHINE SYSTEM 


The simplest form of the swing equation arises from 
consideration of one finite machine, connected through a 
transmission line to an infinite bus. As observed in Example 
3.1, this leads to an equation 


a“ 6 


M—~s~ = P.- P. sin 6= _ £(6) (By L) 
at il m 


Derivation of this equation assumes no damping and constant 
internal voltage as well as constant input. If damping and 
field variation are taken into account, there results the 


equation 
M— y, = f(t, 6, —) (By 2) 


For a multimachine system, the power equation for 


the ke} machine, derived using the assumptions of (7.1),1i1s 
2 


— hike “pif 
M, ——- = P = E.nSey cos (86 - Sy + 6) 
k ar? 11k 4=1 Jos Ss k 5 
= F(o,, So, ‘ She ‘ 5) (Bas) 
Here E is the voltage behind transient reactance of the 


Aha machine; 


is the self- or transfer-admittance element of 
the network facing the machine, and 


n is the total number of machines on the system. 
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B=2 


The general equation for the Sat machine, taking into account 


both damping and field decrement, is of the form 


a*s, 
My = f(t,6 


1’ w) (B.4) 


2 cee 


The application of Runge-Kutta to a system of n equations re- 


presenting an n-machine system is as follows, for the eu 
machine: 
First approximation: 8 (x) S 5 (le) 

YK) > evi: oe AWe (Ky } at t = a 
Second approximation: 

Scene oye ok): fers (chs 

Vk) = VOT) + 2(yKy (x)? : at t = ty + ah 
Third approximation: 

eee oie Pode o (ky FP itiik ? “1° 2Kk) 

= ° = + 

V (k) YG (k) + 2(¥1Ky (x) + €1ko (xy) i hog he ty ah 

Fourth approximation: 
= 5 + 
eel othe so Oleh Pol (ky 22th) ) P23. 0k) 
eae n Ce eo ieee 22 Ck) ae eaa () ’ 
rat t.-= ie “f a5h 
where 
he as 
Ratk) - 9m (a9) 
ue op re 


gnvooos odmi poids ,snt 


(h.a) (6 ao > oraalttt r : ‘a 


134 ort sot \ewolfot en. ab 


ie (a) 2% (apo A 64 ( ae it 
lo + oF = 3 38 Filggy tvs * * igh? oa 7 
‘ ose! i: ree ae 


oltsmixoxqgs | tial 
| "ar "; | 7 hau) * 4 : 
cause? * wate 9 — ie 


foe 
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4 

es Toot) mae Oe): & ae Yom*m(k) ! 
4 

Y2¢K) = Yoo) * 2) Tamme) 


The coefficients involved in the calculation of the oka ck) S 


depend on the particular Runge-Kutta variation used. 
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APPENDIX C 


SOME COMPUTER RESULTS 


KUTIA REPEATEO--VARIOUS H~~EQNs (3,21) 


vee sEXAMPLE NOw 


eeeeVALUES OF DELTA ON THIS PAGE 


TIME 


0,000 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150 
0.175 
0.200 
0.225 
0.250 
0.275 
0.300 
0.325 
0.350 
0.375 
0.400 
0.425 
0.450 
0.475 
0.500 
0.525 
0.550 
0.575 
0.600 
0.625 
0.650 
0.675 
0.700 
0,725 
0.750 
0.775 
0.800 
0.825 
0.850 
0.875 
0.900 
0,925 
0.950 
0.975 
1.000 


KUTTA REPEATED--VaARIOUS II--EON, 


H=0,200 


18.105 


50.544 


98.317 


130.847 


109,877 


26.636 


eee cEXAMPLE NO. 


O. 100 


18.105 


27.586 


50.830 


17.368 


100.093 


118.396 


136.016 


145.336 


139.269 


110.657 


40.864 


0.050 


18.105 
20.569 
27.592 
38.173 
50.923 
64.420 
17.527 
69.560 
100.294 
109.868 
116.672 
127.285 
136.461 
144.557 
146.423 
145.891 
141.805 
132.990 
1L7.176 
GL.3l4 


53.635 


eee e VALUES UF UMEGA ON THIS PAGE 


TIME 


0.c00 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150 
0.175 
0.200 
0.225 
0.250 
0.275 
0.300 
0.325 
0.350 
0.375 
0.400 
0.425 
0.450 
0.475 
0.500 
0.525 
0.550 
0.575 


0.600. 


0.625 
0.650 
0.675 
0.700 
0.725 
0.750 
0.775 
0.800 
0.825 
0.850 
0.875 
0.900 
0.925 
0.950 
0.975 
1.000 


420.200 


0.000 


32.337 


23.929 


204904 _ 


=43.098 


~69.208 


0.100 


0.000 


22.479 


33.291 


314543 


25.199 


21.379 


24.135 


1.992 


-16.377 


-57.961 


-LLL.519 


0.050 


0.000 
12.159 
22.533 
29.783 
33.367 
33.642 
31.608 
28.468 
25.266 
22.767 
21.498 
21.883 
24.409 
11.902 

2.788 
-5.501 

=15.416 
=29.617 
=50.795 
=79.443 


-107.323 


1, CLEARING TIME 0.60 SECONDS. 


0.025 


18.105 
16.727 
20.569 
23.561 
27.594 
32.522 
38.177 
44.375 
50.930 
57.667 
64.429 
1.086 
17.537 
83.712 
A9.5T1L 
95.099 
100.306 
105.219 
109,681 
114.348 
116.687 
122.976 
177.303 
Lil. 768 
146.484 
140.661 
143.593 
195.486 
146.474 
146.631 
145,968 
144.437 
141.925 
136.257 
133.184 
126.387 
117.486 
196.072 
91.706 
74.448 
54.260 


(3.21) 


Ly CLEARING TIME 0.60 SECUNDS. 


0.025 


0.000 
6.197 
12.160 
17.671 
22.537 
26.608 
29.787 
32.034 
33.370 
33.869 
33.645 
32.841 
JLe6lk 
30.107 
28.470 
26.822 
25.269 
23.895 
22.770 
21.956 
21.503 
21.463 
21.891 
22.851 
242422 
17.573 
11.922 
T.1LL7 
2.830 
-1.249 
-54420 
-9.988 
-15.273 
-21.618 
-29.362 
-38.909 
50.449 
63.996 
-79.032 
-F46191 
-107.054 


0.0125 


18.105 
18.727 
20.569 
23.561 
27.594 
32.523 
38.178 
44.375 
50.931 
57.668 
64.430 
71.087 
77.538 
63.713 
89.572 
95.100 
100.307 
105.220 
109.881 
114.348 
118.688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.48K 
146.477 
146.634 
145.972 
144.441 
141.931 
138.264 
133.193 
126.398 
117.499 
106.089 
91.806 
74.6473 
54.287 


0.0125 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32,035 
33.371 
33.869 
33.645 
32.841 
31.611 
30.107 
28.470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.922 
7.118 
2.831 
-1.246 
5.416 
-9.983 
-15.266 
-21.609 
-29.371 
-38.896 
50.434 
=63.979 
-79.014 
“94,175 
-107.043 


0.00625 


16.105 
16.727 
20.569 
23.56) 
27.594 
32.523 
38.178 
44.375 
50,931 
57.668 
644430 
71.087 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109,881 
114.348 
118.688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.488 
146.477 
146.634 
145.972 
144.441 
141.931 
136.264 
133.193 
126.398 
117.499 
106.089 
91.806 
74.472 
54.287 


0.00625 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.841 
31.611 
30.107 
26.470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
244423 
17.574 
11.922 
7.118 
2.831 
1.246 
-5.416 
-9.983 
-15.266 
-21.610 
-29.372 
38.897 
-50.434 
-63.979 
-79.015 
-94.176 
-107.043 


2003125 


18.105 
18.727 
20.569 
23.561 
27.594 
32.523 
38.178 
44.375 
50.931 
57.668 
64.430 
71.087 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109.881 
114.348 
116.6848 
122.977 
127.304 
131.769 
136.486 
140,662 
143.595 
145.487 
146.476 
146.634 
145.971 
144.441 
141.930 
138.263 
133.191 
126.396 
117.497 
106.087 
91.803 
74.469 
54.283 


2003125 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
34.045 
32.041 
3Leoll 
30.107 
28.470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.922 
72118 
2.831 
-1.247 
5.417 
-9.983 
“15.267 
“21.611 
-29.373 
-38.899 
-50.436 
-63.982 
-79.017 
-94.178 
-107.045 
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GILL REPEATED~-VARIUUS H=-EQN. (3.23) 
eeaeEXAMPLE NOw Ly CLEARING TIME 0460 SECONDS. 


aeee VALUES OF DELTA ON THIS PAGE 


TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 +003125 


0.000 18.105 18.105 18.105 16.105 18.105 16.105 16.105 


0.025 18.727 16.727 18.727 16.727 
0.050 20.568 20.569 20.569 20.569 20.569 
0.075 , 23.561 23.561 23.561 23.561 
0.100 27.582 27.592 274994 27.594 276594 27.994 
O.l25 32.522 324523 326523 32.523 
0.150 SBeL?2 38.177 38.178 38.178 36.178 
O.1TS 44.375 44.375 44.375 44.375 
0.200 50-251 50.419 50.922 50.930 50.931 50.931 50.931 
0,225 . 57,667 57.668 57.668 57.668 
_ 04250 64.419 64.429 64.430 64.430 64.430 
0.275 71.086 71.087 71.087 71.087 
0.300 77.348 77.526 77.537 77.538 77.538 77.538 
0.325 63.712 83.713 83.713 63.733 
0.350 89.559 89.571 89.572 89.572 69.572 
0.375 95.099 95.100 95.100 95.100 
0.400 97.598 100.067 100.293 100.306 100.307 100.307 100.307 
0.425 105.219 105.220 105.220 105.220 
04450 109.866 109.881 109.881 109.881 109.981 
0.475 114.348 114.348 114.348 114.348 
0.500 118A, 360 118.670 110.687 LL8.688 116.688 L16.668 
0.525 122.976 122.977 122.977 122.977 
0.550 127.283 127.503 127.304 127.304 127.304 
O.575 31.7698 131.769 131.749 131.769 
0.600 129,362 135.999 136.459 136.484 136.496 136.486 136.486 
0.625 140.661 140.662 140.662 140.662 
0.650 143.553 144.593 143.595 143.595 143.555 
0.675 145.405 145.488 145.468 145.487 
0.700 145.187 146.416 146.474 146.477 146.477 146.476 
0.725 146.631 146.634 146.634 146.034 
0.750 145.879 145.9608 145.972 145.972 145.971 
0.77 144.430 144.441 144.441 144.441 
0.800 100.074 139.820 ‘oLl.783 141.924 141.931 141.931 141.930 
O.A25 138.256 138.264 138.264 136.263 
0.850 132.953 133.182 133.192 133.192 133.191 
0.875 126.364 126.398 126.398 126.396 
0.900 109.404 Lh7.ti3 117.482 117.499 LL? +99 LL7.497 
0.925 196.068 106.089 106.089 106.087 
U.950 91.215 91.780 91.806 91.806 91.003 
0.975 74.441 74.472 74.472 742469 
1.c0o0 -44.415 38.404 945.498 54.252 54.247 54.287 54.283 


GILL REPEATEN=-VARIUUS H=-ENNe (3423) 


eee eEXAMPLE AG. Ly CLEARING FIME 0.60 SFCONDS. 


eeee VALUES OF OMEGA ON THIS PAGE 


TIME He0.200 0.100 0.050 0.025 0.0125 0.00625 2003125 
0.000 0.0C0 2000 0.900 0.000 0.000 0.000 0.000 
0.025 6.197 6.197 6.197 6.197 
0.050 12.158 12.160 12.161 12.16l 12.161 
0.075 17.671 17.671 17.671 17.671 
0.100 22.470 220533 220837 226537 22.837 22.537 
0.125 : 266608 26.608 26.608 26.608 
0.150 200783 296 TAE = 24a KT = 294787 = 294 787 
0.175 Svaba4 sd2s0a5 .32,035 4.392.035 
Beene = 310970, s3s279° “33.986 GoaeavO Peas Ti TSa7t | e23parl 
06225 34.669 336869 336869 33.469 
0.250 4.641 33.645 43.645 33.685 33.645 
0.275 32.841 32.841 32.841 324841 
0.300 31.535 41.408 31.611 Pest iosIGGll. aga lsbil 
0.325 30.100 30.107 30.107 304107 
04390 2be4o? PH 4O9 = 2BeAZO = 284470 28470 
0.375 2busee 26.822 26.822 26.822 
0.400 236704 = 25411 25.266 250269 256269 25.269 25.269 
Oehed 230094 230895 =. 230895 = -2348 95 
0.450 Oe 76? 5 eaTYO seeedtl Sezer See. 77l 
0647S 210956 216956 216956 = 214.956 
0.500 2te36t 2Led97 «= Be SO3. «= 210503 = 216503 21.2503 
0.525 ZLe463 2146463 210463 210 43 
04550 21.882 210891 210892 21.892 21.892 
hah Oe Tee ys Het tt 
24,095 246407 444 . . . 
Sond oy ; LteS73LTeS74—ATABTA TASS. 
04640 11.897 11.921 Migaze CoWcSee Gal lsyee 
0.675 Tells 7.118 7.118 7.118 
0.700 1.782 2.778 2.829 2.831 2.831 2.831 
0+725 10250 10286 3 = = ho BbG |= Le 247 
0.750 -5.519 5.421 50416 = 50416 = 50 417 
aoe Sucas ciracepe otaebs erucus? 
~55. - 2983 -15.446 “15. “15.6 -15. “1d. 
eet ae “7 m210620  =210610 21.610 =21.611 
0.850 290667 = 2963H5. = 296372 296372 «= 296373 
Te Stee eliiees ee 
= -5O.H72 -50.453 e . -_ . = . 
an aos230 =64,00L 63.979 63.979 63.982 
0.950 -19.543 =79,038 -79,015 =79.015  =79,017 
0.975 -94.197 940176 940176 =94.178 


1.000 =500567 -L1?26447 =1076402 =1U7405¥ =107.043 =107.043 =107.045 
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STRACHEY REPEATED=-VARIGUS H-~EQN. (3.25) 
sascEXAMPLE NO« 1, CLEARING TIME 0460 SECONDS. 


aseeaVALUES OF DELTA ON THIS PAGE 
TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 +003125 


0.000 18.105 18.105 18.105 18,105 18,105 18.105 18.105 


0.025 16.727 18.727 18.727 18.727 
0.050 20.569 20.569 20.469 20.569 20.569 
0.075 23.561 23.561 23.561 23.561 
0.100 27.595 274592 27.594 27.594 27.594 27.594 
0.125 32.523 32.523 32.523 32.523 
0,150 38.174 38.177 38.178 38,178 38.178 
0.175 44.375 44.375 44.375 44.375 
0.200 51.449 50.862 50.924 50.930 50.931 50.931 50,931 
0.225 57.667 57.668 57.668 57.668 
0.250 64.422 64.429 64.430 64-430 644430 
0.275 : 71.086 71.087 71.087 71.087 
02300 T1419 17.529 17,537 77,538 77,538 77.538 
0.325 83.712 83.713 83.713 83,713 
0.350 89.562 89.571 89,572 89,572 89.572 
0.375 A 95.100 95.100 95.100 95.100 
06400 100.541 100.159 100.297 100.306 100.307 100.307 100.307 
0.425 105.219 105.220 105.220 105.220 
0.450 109.871 109,881 109.881 109,861 109,881 
0.475 1146348 = 1146348 = 40348 = 1144348 
0.500 116.487 116.676 116.687 116.688 118,688 116.688 
0.525 122,976 122,977 122.977 122.977 
0.550 127.290 1276304 127.304 127.304 127.304 
0.575 131.768 131.769 131.769 131.769 
0.600 135.552 1366163 1366468 136.4465 1366486 13664866 136.486 
0.625 140.661 140.662 140.662 140.662 
0.650 143.567 1434594 143.595 143.595 143,595 
0.675 145.486 145.488 145.488 145.487 
0.700 1456716 = 1466440 = 1464475 = 1464477 = 1466477 = 1464476 
0.725 146.632 «146.634 146.634 146.634 
0,750 145.922 145.970 145.972 145.972 145.971 
0.775 1446439 144.441 144.441 1440441 
0.800 138.992 1406413 141.858 141.928 141.931 141.931 141.930 
0.825 138.261 138.264 138.264 138.263 
0.850 133,082 133,189 1334193 133.193 133.191 
0.875 126.393 126.398 =6126.398 =: 126. 3396 
0.900 113,661 1176329 11746494 117-500 117.499 117.497 
0.925 106.083 106.090 106.089 106.087 
0.950 91.555 91.799 91.807 91.806 91.803 
0.975 74.464 74.474 74.472 74.469 


1,000 796133 47.307 53.969 54.279 54.289 54.287 54.283 


STRACHEY REPEATED=-VARIUUS H=-EUN. (3.25) 
see sEXAMPLE NOs Ly CLEARING TIMF 0.60 SECONDS. 


eeeeVALUES OF UMEGA ON THIS PAGE 


TIME H=0.200 0.100 0.050 0.025 C.0125 0.00625. 003025 
0.C00 0.000 0.000 0.000 0.000 0.000 0.000 06.000 
0.025 6.197 6.197 6.197 6.197 
0.050 2.189 120160 2a OL «= ARB = AZ LO 
0.075 Vibert Sf S7ch7h | Sh %eby1) Ch Yeael 
0.100 22.501 22.534 22.537 22.537 22.537 22.537 
0.125 26.608 26.608 26.608 26.608 
0.150 29.785 29.787 29.787 29.787 29.787 
0.175 32.034 32.035 32.035 32.035 
0,200 «330468 «= 330329 «330908 «= 3.36370 «= 33437L «= -33637L «= 330 BT 
0.225 33.869 33.869 33,869 33,869 
0.250 33.643 33,645 336685 33.645 33.648 
0.275 320841 3208041 326841 324841 
0.300 Sivhha 7H03Gbo.. "OUsAak *SLGORL BLeRha [dk eil 
0.325 30.107 30,107 30.107 30,107 
0.350 28.468 286.470 24.470 28.470 26.470 
0.375 26.822 aeitce Se cae HEE 
4? «-BSe22k 25.267 = 25269 5.26 : : 
Het eo 23.895 23.895 23.895 23.895 
0.450 220769 o22eTTL «—«-so22e TTL —«-s220TTL «22 TTL 
0.475 21.956 214956 214956 214.986 
0.500 210425 geS00. 216803. 214503214503 21.503 
0.525 21.463 21.463 21.463 21.463 
0.550 2ieane = - 216691 216892 = 216892 21.892 
HH SAE tH OR PT 
‘ 7 24.413 . . . . 
Til RAsieR aE - 1te574 170574) ATASTH TAS TS 
04650 11e916 0 LLe922.sALL 922s GAARA 
a ite Ge ORM Bat 
2.613 . . . . 
Sires agers =1c247? «= 1628601024600 10287 
0.750 =5.458 9 =5.418 50416 9 = Behl = = 50 417 
niet Teizes “albctes 5.860 -ipcRe? 
= = =15.344 -15.269 = . = . = . 
ping ta aah =21.613 21.609 216609 =21.611 
0.850 29.496 29,375 =294371 «29.372 © = 294373 
st a oe she ae 
s = =50.438 ben . = . = . 
ae Perel eeee ee eggs =63,978 83.079 03.982 
0.950 -79.197 =79,018 =79.013 =79.015 =79,017 


0.975 ~O4,177 =94.174 <=94.176 =94.178 
1.000 63.239 =108.802 <107.129 =107,043 =107.042 ~107.043 =107.045 


BOULTON REPEATED-=VARIOUS H=-EON. (3.27) 
seaeEXAMPLE NOw Ly CLEARING TIME 0.60 SECONDS, 


aososVALUES OF DELTA ON THIS PAGE 
TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 .003125 


0.000_ 18.105 18.105 18.105 18.105 18,105 18.105 16,105 


0.025 18.727 18.727 18.727 16.727 
0.050 20.569 20.569 20.569 20.569 20.569 
0.075 23.561 23.561 23.561 23.561 
0.100 27.589 27.592 27.594 27.594 27.594 27.594 
0.125 32.522 32.523 32.523 32.523 
0.150 38.173 38,177 38.178 36.178 38.178 
0.175 44,375 44.375 44.375 44.375 
0.200 50.817 50.843 50.923 50.930 50.931 50.931 50.931 
0.225 ‘ 57.667 57.668 57.668 57.668 
0.250 64.421 644429 644430 64.430 64.430 
0.275 71.086 71,087 71.087 71.087 
0.300 717.383 77.527 17,537 77.538 77.538 77.538 
0.325 63.712 63.713 83.713 83.713 
0.350 89.560 89.571 89.572 89.572 89.572 
0.375 95.099 95.100 95.100 95.100 
0.400 98,042 100.100 100.294 100.306 100.307 100.307 100.307 
0.425 105.219 105.220 105.220 105.220 
04450 109.867 109.881 109,881 109.8al 109,881 
0.475 114.348 114.348 114.348 114.348 
0.500 118.393 118.671 118.687 118.668 118.688 116.688 
0.525 122.976 122.977 122.977 122.977 
0.550 127.284 127.303 127.304 127.304 127.304 
0.575 131.768 131.769 131.769 131.769 
0.600 130,003 136.004 136.460 136.484 136.486 136.486 136.486 
0.625 140.661 140.662 140.662 140.662 
0.650 143.553 143.593 143.595 143.595 143.595 
0.675 145.485 145.468 145.488 145.487 
0.700 145.267 146.416 146.474 146.477 146.477 146.476 
0.725 146.630 146.634 146.634 146.634 
0.750 145.860 145.967 145.972 145.972 145.971 
0.775 144.436 144.441 144.441 144.441 
0.800 105.537 139.115 141.785 141.924 141.931 141.931 141.930 
0.825 136.255 138.264 138.264 136.263 
0.850 132.957 134.162 133.192 133.192 133.191 
0.875 126.384 126.398 126.398 126.396 
0.900 110.296 117.124 117.482 117.499 117.499 117.497 
0.925 106.068 106.089 106.089 106.087 
0.950 91.238 91.781 91.606 91.606 91.803 
0.975 74.442 74.472 74.472 74.469 


1.c00 -19.877 406474 53.540 54.253 54.287 54.267 54.283 


BOULTCN REPEATEO--VARICUS H--EQN. (3.27) 
eee sEXAMPLE NCs Ly CLEARING TIMF 0.60 SECUNDS. 


eeeeVALUES OF OMEGA ON THIS PAGE 


TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 .003125 
0.000 0.000 0.000 0.000 0.000 0-000 0.000 0.000 
0.025 6.197 6.197 6.197 66197" 
0.050 12.156 12.160 12.161 12.161 L2.1él 
0.075 17.671 17.671 17.671 17.071 
0.100 22.470 22.533 224531. . 226537 220837, 226537 
0.125 26.608 264608 26.608 264608 
0.150 29.783 29.787 29.787 29.787 29.787 
0.175 32.034 32.035 32.035 324035 
0.200 31.887 334287 33.367 33.370 334371 33.371 33.371 
0.225 33.869 33.869 33.869 334869 
0.250 33.642 33,645 33.645 33.645 33.645 
0.275 32.941 32.841 32.841 32.841 
0.300 314544 314608 31.611 31.611 31.611 31.611 
0.325 30,107 30.107 30.107 304107 
0.350 28.467 28.469 28.470 284470 28.470 
0.375 26.822 264822 26.822 26.822 
0.400 23.794 25.195 254266 25.269 254269 254269 25.269 
0.425 23.895 234895 23.895 234895 
04450. 22.767 P2avi0in. 224771 22.771 22.771 
0.475 21.956 214956 21.956 214956 
0.500 21.371 21.497 21.503 214503 21.503 21.503 
04525 21.463 21.463 21.463 2146463 
0.550 21.882 21.891 214892 21.892 21.892 
5 2 ou eter (OMe eke 
OVID S 240eOT se 2408 he ; : 
te eg es B hath, ‘lishtea Like. “Losts 
0.650 11.898 11.921 11.922 11.922 11.922 
0.675 7.116 7.118 7.118 7118 
0.700 1.944 2.780 2.829 2.831 2.831 2.831 
0.725 14290, <10246. <2e246, hed? 
0.750 =5.517 -5e42L -50416 54416-52417 
“0.775. -9.989 -9.983 -9.983 -9.983 


- - 61 -154444 -15.275 -15.266 -15.266 -15.267 
oiets rahi se : -21.620 -21.610 -21.610 21.611 


i = 290663 9 --294385 «9 -294372 «294372 292373 
0.878 =38.913 -38.897 marth Hits 
0. =58.545 -506869 -506454 -50+6434 -50.434 -50+436 
PL apenas? = 64.002 -63.979 63.980 63.982 
0.950 =719.540 -79.038 -79.015 -79.015 <-79.017 
0.975 =94.198 94.176 -94.176 -94,178 


1.000 44.419 -111.726 -107.383 -1074059 -107.043 ~-107.043 107.045 


= se We al 
- se 
= ey ; 
a _ 
i —~, 
: f se  ¢ 
& : a 
> - ie . 


= s 
= >. af a 
; 7 a : / 
i ‘ 
¢ - a = 7 cs 
= ° - 
§ a + - | 4, 
a = a E ee Vee a e ba i 
bat te 7 : ? : oe Ney : 
: : >. s 
a. a (fi 7 i. pa , - 
- 7 = - ~~ a P 7 
: i an ae 
- Pi Fe Dam ee 
' - 7 
y , r + > ball 5 3 
=e* 7s : mich 14 ee : 
7 a. | a —- » = » 
, — ; < % xy om Me t-e> te ~ “x D 
7 Z z=! me Pet ing ae J « Arad 
2 ’ : 
i ache ’ i. ne oes om 
— 3 Sh a, ny ie ai 
of be ~~ a ek ‘ee ; en 
ball £Z — de S | ese a si sw aa - 


Teeter tee et teense ees s 


MERSON REPEATEOD=-VARIOUS H--EQN. (3.34) 
s¢eeEXAMPLE NO. Ly CLEARING TIME 0.60 SECONDS. 


-seeeVALUES OF DELTA ON THIS PAGE 


TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 +003125 


0.C00 18.105 16.105 16.105 14.105 16.105 18.105 16.105 


0,025 186.727 18.727 16.727 18.727 
0.050 20.568 20.569 20.569 20.569 20.569 
0.075 23.561 23.561 23.561 23.561 
0.100 27.582 27.593 27.594 27.594 27.594 27.594 
0.125 32.523 32.523 32.523 32.523 
0.150_ 38.176 38.178 38.178 38.178 36.178 
0.175 44.375 44.375 44.375 44.375 
0.200 50.196 50.694 50.929 50.931 50.931 50.931 50.931 
0.225 57.668 57.668 57.668 57.668 
0.250 644428 64.430 64.430 64.430 64.430 
0.275 71,087 71.087 71.087 71.087 
0.300 77.500 77.536 77.538 77.538 77.538 77.538 
0.325 63.713 83.713 83.713 83.713 
0.350 89.570 89.572 89.572 89.572 89.572 
0.375 95.100 95.100 95.100 95.100 
0.400 99.550 100.281 100.306 100.307 100.307 100.307 100.307 
0.425 105.220 105.220 105.220 105.220 
0.450 109.861 109.881 109.881 109.861 109.881 
0.475 114.348 114.348 114.348 114.348 
0.500 118.663 118.688 118.688 118.688 116.688 118.688 
0.525 122.977 122.977 122.977 122.977 
0.550 127.304 127.304 127.304 127.304 127.304 
0.575 131.769 131.769 131.769 131.769 
0.600 135.127 136.458 136.486 136.486 136.486 136.486 136.486 
0.625 140.662 140.662 140.662 140.662 
0.650 143.595 143.595 143.595 143,595 143.595 
0.675 145.488 145.488 145.488 145.487 
0.700 146.439 146.476 146.477 146.477 146.477 146.416 
0.725 146.635 146.634 146.634 146.634 
0,750 145.975 145.972 145.972 145.972 145.971 
0.775 144.442 144.441 144.441 144.441 
0.800 134.919 141.840 141.937 141.932 141.931 141.931 141.930 
0.825 138.266 138.264 136.264 138.263 
0.850 133.204 133.195 133.193 133.193 133.191 
0.875 126.400 126.398 126.398 126.396 
0.900 117.260 117,519 117.502 117.500 117.499 117.497 
0.925 106.093 106.090 106.089 106.087 
0.950 91.832 91.811 91.808 91.806 91.803 
0.975 74.478 74.474 74.472 7142469 
1.000 11.412 53.581 54.319 54.293 54.289 54.287 54.283 


MERSON REPEATEO--VARICUS H--EQN. (3434) 


eeecEXAMPLE NOw Le CLEARING TIMF 0.460 SECONDS. 


éeeeVALUES OF UMEGA ON THIS PAGE 


TIME H=0.20U 0,100 0.050 0.025 C.0125 0.00025 .003125 
0.000 0.060 0.000 0.000 0.000 0.000 0.000 0.000 
0.025 6.197 6.197 6.197 6-197 
0.050 12.160 12.160 12.161 12-161 12.161 
0.075 17.671 17-67? 17.672 17.671 
0.100 22.528 224536 22.537 22.537 22.537 22.537 
0.125 26.608 26.608 26.608 26.608 
0.150 29.787 29.787 29.787 29.787 292787 
0.175 32.035 32.035 32.035 32.035 
0.200 33.151 33.369 33.371 33.371 33.371 33.371 33.371 
0.225 33.869 33.869 33.869 33.869 
0.250 33.645 33.645 33.645 33.645 33.645 
0.275 32.841 32.841 32.841 32.841 
0.300 : 31.618 31.611 31.611 31.611 31.611 31.611 
0.325 30.107 30.107 30.107 30.107 
0.350 26.470 26.470 28.470 26.470 28.470 
0.375 | 26.822 26.822 26.822 26.822 
0.400 25.131 254274 25.269 254269 25.269 25.269 25.269 
04425 i 23.895 23.895 23.895 23.895 
0.450 22.771 22.771 22.771 22.771 22.771 
0.475 21.956 21.956 21.956 21.956 
0.500 21.501 21.503 21.503 21.503 21.503 21.503 
0.525 214463 21.463 21.463 216463 
0.550 21.892 21.892 21.892 21.892 21.892 
os203 gg A 
4.408 24.423 246423 4s a ‘ 
a eas poy ae : 17.574 17.574 17.574 17.574 
0.650 11.924 11.923 11.922 11.922 11.922 
0.675 7.118 7.118 7.118 7.118 
0.700 2.782 2.835 2.832 2.832 2.831 2.831 
0.725 -1.246 -1.246 1.246 -1.247 
0.750 =5.411 -5.415 5.416 5.416 5.417 
0.775 -9.982 -9.982 -9.983 -9.983 


. -26.002 -15.377 -15.256 -15.265 -15.266 -15.266 -15.267 
0.825 ; -21.607 -21.609 -21.609 -21.611 


0.850 290353 294369 «290371 «90 290372 © = 290373 
0.675 -38.893 -38.896 bet he 

=50.630 -502404  -50.430  -50.433. -50+434  -50+436 
ears ee raae 63.974 636978 63.979  -63.982 
0.950 -78.975  -79.009 -794013 79.015 -T19.017 
0.975 =-94.170 —-946174 94.176  -94.178 


1.000 -108.254 -106.998 -107.010 -107.039 -107.042 -107.043 -107.045 


MERSON REPEATEO--VARIOUS H=-EUN. (3.34) 
aeeeEXAMPLE NOw Ly CLEARING TIME 0.60 SECONDS. 


eeeeERROR IN CELTA ON THIS PAGE, VALUES*#1000. 


TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 .003125 
0.000 0.000 7200.042 -96.411 -18.716 1.425 -0.083 -0.005 
0.025 -0.000 -0.000 -0.000 -0.000 
0.050 ~0.006 -0.004 -0.000 -0.000 -0.000 
0.075 -0.010 -0.001 -0.000 -0.000 
0.100 -94426 -0.257 -0.018 -0.001 ~0.000 -0.000 
0.125 ~0.029 -~0.002 -0.000 -0.000 
0.150. -0.623 -0.041 -0.003 -0.000 -0.000 
0.175 -0.054 -0.003 -0.000 -0.000 
0.200 -32.050 <-15.089 ~1.032 -0.066 -0.004 -0.000 -0.000 
0.225 -0,.075 -0.005 -0.000 -0.000 
0.250 -1.310 -0.081 -0.005 -0.000 -0.000 
0.275 -0.082 -0.005 -0.000 -0.000 
0.300 ~23.592 -1.344 -0.080 -0.005 -0.000 -0.000 
0.325 -0.075 -0.005 -0.000 -0.000 
0.350 ~1.180 -0.068 -0.004 -0.000 -0.000 
0.375 -0.061 -0.004 -0.000 -0.000 
0.400 -410.015 <-18.077 -0.948 -0.054 -0.003 -0.000 -0.000 
0.425 -0.048 -0.003 -0,000 -0.000 
0.450 -0,748 -U.042 -0.003 -0.000 -0.000 
0.475 ~U.038 -0.002 -0.000 -0.000 
0.500 -11.774 -0.619 -0.035 -0,002 -0.000 ~0.000 
0.525 -0.033 -0.002 -0.000 -0.000 
0.550 -0.569 -0.033 -0.002 -0.000 -0.000 
0.575 -0.033 -0.002 -0.000 -0.000 
0.600 ~-195,456 10.575 -0.601 -0.035 -0.002 -0.000 -0.000 
0.625 -0.088 -0.004 -0.000 -0.000 
0.650 1.752 -0.051 -0.002 ~0.000 -0.000 
0.675 -0.028 -0.00L -0,000 -0.000 
0.700 -41.893 -0.695 -0.014 -0.000 -0.000 -0.000 
0.725 -0.005 -0.000 -0.000 -0.000 
0.750 -O0.141 0.U01 0.000 -0.000 -v.000 
0.775 0.004 -0.000 -0.000 -0.000 
0.800 -1124,845 -4.573 0.156 0.001 -0.001 -0.000 -0.000 
0.825 -0.010 -0,002 -0.000 -0.000 
0.850 0.021 -0.039 -0.004 -0.000 -0.000 
0.875 -0.106 -0.009 -0.001 -0.000 
0.900 3.756 -1.621 -0.219 -0.018 -0.00L -0.000 
0.925 -0.429 -0.034 -0.002 -0.000 
0.950 -7.849 -0.758 -0.057 -U.004 -0.000 
0.975 -1.172 -0.083 -0.005 -0.000 
1.000 -1023.836 -172.U82 -70.314 -1.498 -0.098 -0.006 -0.000 


MERSON REPEATED=-VARIOUS H--EQN. (3.34) 


eeesEXAMPLE AC. Ly CLEARING TIME 0.460 SECUNDS. 


eeeeERROR IN GMEGA ON THIS PAGE, VALUES#1L000. 


TIME H=0.200 0.100 04050 0.025 0.0125 0.00625 .003125 
0,000 0.000 -4904.946 1518.484 66.854 2.000 0.20% 0-015 
0.025 =-0.017 -0,001 -0.000 0.000 
0.050 -0.541 0.027 -0.001 -0.000 -0.000 
0.075 -0.037 -0.002 -0.000 -0.000 
0.100 -17.405 -0.857 -0.048 -0.003 -0.000 -0.000 
0.125 -0.057 -0.004 -0.000 -0.000 
0.150 -1.077 -0.062 -0.004 -0.000 0.000 
0.175 -0.059 -0.003 -0.000 -0.000 
0.200 -523.117 -20.087 -0.415 -0.048 -0.003 -0.000 -0..000 
0.225 -0.03L -0.002 -0.000 -0.000 
0.250 -0.352 -0.012 -0.000 -0.000 -0.000 
0.275 0.007 0.008 0.000 0.000 
0.300 -0.468 04244 0.022 0.002 0.000 0. 000 
0.325 O.031 0.002 0.000 0.000 
0.350 0.554 0.036 0.002 0.000 0.000 
0.375 0.036 0.002 0.000 0.000 
0.400 131.565 9,731 0.568 0.033 0.002 0.000 0.000 
0.425 0.029 0.002 0.000 0.000 
0.450 0.425 0.023 0.001 0.000 0.000 
0.475 0.017 0.001 0.000 0.000 
0.500 5.729 0.238 0.011 0.001 0.000 0.000 
0.525 0.005 0.000 0.000 0.000 
0.550 0.044 -0.001 -0.000 -0.000 0.000 
ets BOIL Mapeaed /caaceaas bonieon 
-0.531 -0.170 -0.014 -0. -0. -U. 
EL —- - 0.223 0.0L 0.001 0.000 
0.650 3.677 0.135 0.006 0.000 0.000 
0.675 0.084 0.004 0.000 0.000 
0.700 12.954 1.748 0.054 0.002 0.000 -0.000 
0.725 0.035 0.001 0.000 0.000 
0.750 0.937 0.018 -0.000 -0.000 0.000 
Gat Soh -E artes eye epee ah et: 
é 0.219 -0.03 -0. -0. =Dé 
beaes ee pa -0.097 -0.009 -0.001 -0.000 
0.850 -1.797 -0.213 -0.017 -0.001 -0.000 
et Sere 
-5 -8.812 -0.789 -0. -0. ey 
metre = oe -1.327 -0.093 -0.006 -0.000 
0.950 25.828 -1.916 -0.127 -0.008 -0.001 
0.975 -2.107  -0.125 -0.007 -0.000 


1.000 -7323-585 -539.591 -29.530 -1.178 -0.044 -0.002 -0.000 


MERSCN AUTCMATIC PROCEDURE, EQN. (3.36) 


EXAMPLE NO. 


LEME s 
SEC 


0.000 
0.025 
0.050 
0.075 
0.100 
Uses 
0.150 
0.175 
0.200 
Ga2eS 
0.2250 
Use 27.9 
0. 300 
0.325 
0.350 
OeoTa 
02400 
0.425 
0.2450 
Ue475 
0.500 
0.525 
0.550 
0.6C0 
0.625 
0.650 
0.675 
0.700 
Gai25 
0.750 
Ae he 
0.800 
0,825 
0.850 
Oe875 
0.9C0 
0.925 
0.950 
0.975. 
1-0CO0 


ly 


VALUES OF 
DELTA 


18.105 
MRR HH 
20.568 
HHHeHHR HHH 
HHHHRURH HH 
HHRHHH HE 
38.159 
RHHHH REE 
MRR MH H Ht Ht 
HHUA HHH 
64.397 
HHHHHHH H 
HHH HH 
HHRHUHH HAR 
89.546 
HHH RHR HHH 
HoH HHH HH 
H 3 HSH OH OH tt 
109.864 
HR RMR HK 
HHH HK HH EK 
RHUL HH 
TZLAZUG 
HH RHUL HH 
136.464 
HHRGH HHH 
HHKRH REAR 
HHH HA KH HH 
146.452 
HHHRHHHRH 
145.930 
HHRKHHHE 
HHHHHHH HK 
st Ht He OH Ht OE OO 
133.093 
3 tt He He Ht St ot 
HRHRHHHKHH 
HHHHAMR HH 
91.487 
HHHHAH HHH 


532945 


CLEARING TIME 0.60 SECONDS. 


VALUES OF 
OMEGA 


0.000 

HR HHH 
12.160 
eH HO HO 
HH et Ht Ht 4 
HHH HH HH Ht 
29.7183 
HHHR HR HH 
HHUMRKH HH HK 
HR HH HH tt 
33.648 
HH Ht 3 At HEX Gt 
HHHRHRHH H 
RH MH tH tt Ht He Ht 
28.415 
St He HE Ht HE HE tt 
MRR HHH HH HE 
Hott He HM HH 
recite 
HHie RH Re F 
HH MH HK H 
RUMH MA BH tH 
21.885 
HUA HR RH HK 
Ze 4 lO) 
Hott He OF Ht HEHE HE 
HHH RAH HH 
HK HK HHH 
2.f98 
HHH HH Hide Hf 
—5.4/73 

st HEHE St HH SEH tt 
HM Ht HH 
HR HH HH 
=~29.413 
HK RRR 
HHH 4S HE OH Oe 
HR HH HH 
-79.074 


HHeHRH HTH FT 


-107.021 


AE 


oa | 
0.-10E OO 


a1 Ee 
Rea 


Ce050 
HHH HHH 

C.100 
Ha 
RUHR EH 
HE 

C.100 
RH HH 
HHH HH 
ee es 

0.100 
HHH HH 
HHH RH 
HH HE 

C.100 
HHH 
HHH RE 
Ht at at Ht 

C.100 
HH HH 
HHH HH 
#4 tH aE at 

C.100 
HE Ht 

C.100 
ee a 
HoH Ht 
ee 

0.050 
ee ss 

C.100 
HHH 
HR HG GH 
HHH HH 

C.100 
HEHEHE 
KH RH 
HHH HHH 

C.050 
HoH HHH 


0-050 


02020 


2aeeee 


00140 


eae ave 
aeenae 
saveas 
eaveor 
ee ee 
061.9 


oneea a. 


eeeeee 


#anenae 


O01 69 

eoneee 
aeoeus 
soenea 
00} .0 
eeaeene 
aeowe i 
HTH 
OO1L.9 

epnerna 
O0f.9 

aniene 
waaene 
vaneys 
O¢0.0 


weeeoe 


1 ge what 


~~ 


ere re 


tetoe woe 
cBe.1S 


[1a.As > ly ae rt: ° : 
WR Hae vee ae 7 
al erp lage Le 
‘ ’ o8) ¥ ; - 
‘ ‘ } ; ie 
7. iy, ee a! 
7 


JOHNSON=WARD METHOD--VARIOUS H--EQN. 


aeeeEXAMPLE NO. 


(5.12) 


CLEARING TIME 0.60 SECONDS. 


aeeeVALUES OF DELTA ON THIS PAGE 


TIME 


0.c00 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150 
0.175 
0.200 
0.225 
0.250 
0.275 
04300 
0.325 
0.350 
04375 
0.400 
0.425 
0.450 
0.475 
0.500 
0.525 
0.550 
0.575 
0.600 
0.625 
0.650 
0.675 
0.700 
0.725 
0.750 
0.775 
0.800 
0.825 
0.850 
0.875 
0.900 
0.925 
0.950 
0.975 
1.000 


JCHNSCN-WARD METHOC--V4R IOUS H--EUNe 


H=0.200 


18.105 


51.315 


102.442 


141.646 


160.694 


433.669 


eee ec EXAMPLE NOW 


0.100 


18.105 


27.601 


50.986 


77.717 


160.622 


LL9O.LT7 


137.332 


148.546 


148.371 


136.156 


97.789 


0.050 


18.105 
20.569 
27.595 
38.182 
50.942 
64.450 
717.567 
89.610 
100.354 
109.938 
114.759 
127.395 
136.607 
143.776 
146.7861k 
146.498 
142.835 
134.723 
120.012 
95.682 


59-540 


0.025 0.0125 


18.105 18.105 
18.727 18.727 
20.569 20.569 
23.561 23.561 
27.594 27.594 
32.523 324523 
38.178 38.178 


44,371 44.376 
50.932 50.931 
517.670 57.668 
64.433 642430 
71.0°9 71.087 


17.5%2 77.538 
83.718 83.714 
89.577 89.572 
95.106 95.101 
100.313 100.308 
105.226 105.220 
109.889 109.882 
114.357 114.349 
118.697 118.689 
122.988 122.978 
127.316 127.306 
131.783 131.771 
136.502 136.488 
140.681 140.665 
143.619 143.598 
145.518 145.492 
146.517 146.482 
146.687 146.641 
146.041 145.981 
144.531 144.453 
142.049 141.946 
138.417 138.284 
133.391 133.218 
126.653 126.430 
117.824 1t7.541 
106.496 106.141 
92.304 91.870 


75.063 74.548 
54.956 54.372 
(5.12) 


CLEARING TIME 0.60 SECONDS. 


wees VALUES OF OMEGA ON THIS PAGE 


TIME 


0.CO00 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150 
0.175 
0.200 
0.225 
0.250 
0.275 
0.300 
0.325 
0.350 
0.375 
0.400 
0.425 
0.450 
0.475 
0.500 
0.525 
0.550 
0.575 
0-600. 
0.625 
0.650 
0.675 
0.700 
0.725 
0.750. 
0.775 
0.800 
0.825 
0.850 
0.875 
0-900 _ 
0.925 
0.950 
0.975 
1.C00 


H=0.200 


0.0C0 


34.618 


26.089 


29.200 


55.578 


235.553 


0.100 


0.000 


22.598 


33.541 


31.779 


254435 


21.808 


252128 


5.949 


-6.360 


—274214 


-73.798 


0.050 


0.000 
12.163 
22.546 
29.805 
33.394 
33.669 
31.633 
28.491 
256291 
22.800 
21.545 
21.955 
24.520 
12.152 
34249 
-4.699 
-14.062 
—27.411 
-47.466 
-75.276 


-104.320 


0.025 0.0125 
0.000 0.000 
6.197 6.197 

12.161 12.161 
174672 b7.671 

22.938 22-537 

26.610 26.609 

296789 29-788 

32.037 32.035 

336374 33371 

33.872 33.869 

33-648 93-645 

32-6344 32.841 

FI eOEY 3h-611 

30.109 3O~LOT 

28.472 28.470 

26.825 26.823 

25.271 25.269 


23.898 23.895 
22.774 22.771 
21.960 21.956 
21.508 21.504 
21.470 21.464 
21.900 21.893 
22.862 22.853 
244435 24.424 
17.594 17.576 
11.952 11.926 


7.158 7.123 
2.885 2-838 
-1.176 -1.237 
-5.324 -5.404 
-9.862 -9.967 
-15.111 -15.246 


-21.410 -21.584 
-29.119 -29.339 
38.582 -38-856 
50.054 50.385 
=63.540 -634923 
-78.544  -78.955 
-93.732  -946119 
-106.721 -107.002 


0.00625 


18.105 
18,727 
20.569 
23.561 
27.594 
32.523 
36.178 
44.375 
50.931 
57.668 
64.430 
71.087 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109,881 
114.349 
118.688 
122.977 
127.304 
131.769 
136.486 
140.663 
143.595 
145.488 
146.477 
146.635 
145.973 
144.443 
141.933 
136.266 
133.196 
126.402 
117.504 
106.095 
91.814 
74.481 
542297 


0.00625 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33-869 
33-645 
32-841 
31.60 
30.107 
28.470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.923 
7.119 
2.832 
-1.245 
—52415 
-9.981 
-15.264 
-21.606 
-29.368 
-38.892 
-50.428 
-63.973 
-79.007 
-94.169 
-107.038 


2003125 


16.105 
186.727 
20.569 
23-561 
27.594 
32.523 
36.178 
44.375 
50.931 
57.668 
64.430 
71.087 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109.681 
114.348 
118,688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.488 
146.476 
146.634 
145.972 
144.441 
141.930 
138.263 
133.192 
126.396 
117.498 
106.087 
91.804 
74.470 
54.284 


2003125 


0.000 
6.197 
12.161 
L7.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.841 
31.61 
30.107 
28.470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.8692 
22.852 
24.423 
17.574 
11.922 
7.118 
2.831 
-1.247 
-5-417 
-9.983 
-15.267 
-21.610 
-29.373 
-38.898 
-50.436 
-63.981 
-79.017 
-94.178 
-107.044 


ba ai Feus ii ee 


a at 


COLLATZ PROCEDURE--VARIOUS H=-EQN. (5.17) 
oeesEXAMPLE NOs Ly CLEARING TIME 0.60 SECONDS. 


seeeVALUES OF DELTA ON THIS PAGE 
TIME H=0,200 0.100 0.050 0.025 0.0125 0.00625 .003125 


0.000 18.105 18.105 18.105 18.105 16.105 18,105 18.105 


0.025 18.727 18.727 18.727 18.727 
0.050 20.568 20.569 20.569 20.569 20.569 
0.075 23.561 23.561 23.561 23.561 
0.100 27.582 27.593 27.594 276594 27.594 27.594 
0.125 32.523 32.523 32.523 32.523 
(0.150 38.175 38,178 36.178 38.478 38,178 
0.175 44.315 44.375 44.375 44.375 
0.200 50.196 50.881 50.927 50.930 50.931 50.931 50.931 
0.225 57.667 57.668 57.668 57.668 
0.256 644426 644430 64.430 64.430 64.430 
0.275 71.087 71.087 71.087 71.087 
0.300 (17.482 174534 77,538 77.538 77.538 77.538 
0.325 83.713 83.713 83.713 63.713 
0.350 89.568 . 89.572 89.572 89.572 89.572 
0,375 ; 95.100 95.100 95.100 95.100 
0,400 100.261 100.256 100.304 100.307 100.307 100.307 100.307 
0.425 105.219 105.220 105.220 105.220 
02450 109.879 109,681 109.881 109.681 109,881 
0.475 114.348 114.348 114.348 114.348 
0.500 118.628 118,685 118.688 118.688 118.688 118.688 
0.525 122.977  ld2.977) sbe2097d. 5 L22-977 
0.550 127.301 127.304 127-304 127.304 127.304 
0.575 © 131.769 131.769 1316769 131.769 
0.600 1364240 136.386 136.481 136.406 136.486 136.486 136.446 
0.625 140.662 140.662 140.662 140.662 
0.650 143.586 143.595 143.595 143.595 143.595 
0.675 145.4687 145.488 145.488 145.487 
0.700 146.143 1466461 = 1466476 = 1466477 = 1460477 = 1466476 
0.725 146.633 146.634 146,634 146.634 
0.750 145.946 145.971 045.972 245.972 145.971 
0.775 144.440 144.441 1440441 1444441 
0.6800 137.499 141.089 141.993 141,930 141.931 141.931 141.930 
0.825 138.263 138.264 138.264 138.263 
0.850 133.130 gobs 4.19) eels. 199 eel35e00s 10139. 191 
0.875 126.395 126.398 126,398 126.396 
0.900 115.254 117.399 117.496 117.500 117.499 117.497 
0.925 106.085 106,090 106.089 106.087 
0.950 91.648 91.801 91.807 91.806 91.803 
0.975 74.466 74.474 74.472 74.469 
1.c00 30.438 49.436 54.073 54.280 54.289 54.287 54.283 


COLLAT? PROCEDURE--VARTUUS H--EON. 15.17) 


eee sEXAMPLE NOe Ly CLEARING TIME 0,60 SECONDS. 


eee VALUES OF UMEGA GN THIS PAGE 


TIME  H20.200 04100 0.050 0.025 0.0125 0.00625 .003125 
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
0.025 6.197 6.197 6.197 6.197 
0.050 12.161 12.161 12.161 pratal 12.161 
0.075 17.671 L7eoal 17.671 17.67h 
0.100 Bo Sh 7RES: 22-8308u) 220587 Pas 220G3T1a4, 226537 ae 226537 
0.125 26.608 26.608 26.608 26.608 
0.150 290789. 29.787 = 296787 «= 296787 «29 787 
0.175 32.035 32.035 320035 32.035 
D200 dabble” 33cdbaees 996372245. 33.97igae 33497 330d7L cy 330371 
0.225 33.069 33.869 33.869 33.869 
0.250 33.646 «334645 330645 354645 33.245 
0.275 326841, 32684] 326841. 32.841 
0.300 4) pits 21 ced dete Slsbkdedn Bich] ign JlaBl Lo, 32.611 
0.325 30.107 30.107 30.107 30.107 
0.350 28.470 28.470 28.470 286470 28.470 
0.375 260B22... 260822. 26822 26.622 
0.400 256281 254264 250269254269 = 254269 = 254269 254269 
0.425 ; 23.895 234.895 23.895 23.895 
0.450 Jogi ieee 2rattiing 2200 lca.. 2207 tine 22670) 
0.475 Dlcohtiae 2lsvShie 2legsG ee 21.956 
0.500 210492216503. 216503 216503 = 214503 216503 
0.525 Zigeehn Cickhaee Ziceeae 2lss63 
0.550 210892 216892212892 «210892 ~—216892 
0.575 22.852 226852 226852 22.852 
0.600 242649 244389 246422 244423 240423 240423 244423 
0.625 << 17.574 17.574 17.574 17.574 
0.650 11.918 116922 116922 LLe922 e922 
0.675 7.118 7.118 7.118 7.118 
0.700 2.527 2.818 2.831 2.832 2.831 2.831 
0.725 a{,oetee  =lhhGan <“1etheee =1s 247 
0.750 =5.445 Sahl? 50416-54416 = 50417 
aise! Targshal ~ialggeee -15cashe) 150267 
-2 = 4 -15.318 -15.26 = . i . - . 
ht P4 ee ee =2le612  -210609 -21.609  -21.611 
0.850 =29.455 -290374 -2903TL «290372. 290373 
0.875 -38.900 -36.896 -38.897  -38.899 
0.900 -53.038 -50.551 -506437 504433 504434 — - 50.436 
0.925 =63.983 63.978 -63.979  -63.982 
0.950 ~19.146 -79018 -792013 -79.015 -79.017 
0.975 =94.179  -940175 -94.176 94.178 


1.000 -100.361 -108.824 -107.144 -107.046 -107.042 -107.043 -107.045 


SECOND ORDER METHOC=-VARIOUS H--EQN. (5.25) 
o4eeEXAMPLE NOs 1y CLEARING TIME 0.460 SECONDS. 


seeeVALUES OF DELTA ON THIS PAGE 
TIME H=0.200 0.100 0.050 0.025 


0.000 18.105 18.105 18.105 18.105 


0.025 18.727 
0.050 20.569 20.569 
0.075 23.561 
0.100 27.594 27.594 27.594 
0.125 32.523 
04150 38.178 38.178 
0.175 44,375 
0.200 50.945 50.932 50.931 50.931 
0.225 57.668 
02250 64.430 64.430 
0.275 71.087 
0.300 17.549 77.538 717.538 
0.325 63.713 
0.350 89.5973 69.572 
0.375 95.100 
0.400 101,123 100.333 100.308 100.307 
0.425 105.220 
0-450 109,883 109,881 
0.475 114.349 
0.500 118.729 118.690 118.688 
0.525 122.977 
0.550 127.307 127.304 
0.575 131.769 
9.600 138.571 136.552 136.489 136.486 
0.625 140.663 
0.650 143.600 143,595 
0.675 145.4846 
0.700 146.606 146.494 146.477 
0.725 146.635 
0.750 145.983 145.973 
0.775 144.442 
0.800 155.938 142.262 141.9468 141.932 
0,825 : 138.266 
0.850 133.221 133.195 
0.875 126.401 
0.900 118,365 117.544 117.503 
0.925 106.094 
0.950 91.869 91.812 
0.975 74.479 
1.000 143.691 55.841 54.367 542294 


SECOND GRNER METFUL--VARTOUS H--EQNs (5625) 


eeeeEXAMPLE ACo Ly CLEARING TIME 0.69 SECUNDS. 


eeee VALUES OCF OMEGA UN THES PAGE 


TIME H=0.200 0.100 0.050 0.025 

0.000 0.000 6.000 0.00 0.000 
0.025 6.197 
0.050 i2.t6l L2.161 
0.075 17.071 
0.100 22.557 22.538 22.537 
0.125 266608 
0.150 24.789 29.787 
O.17y 32.03% 
0.200 33.933 33.393 33.372 33.371 
0.225 33.869 
0.250 33.646 33.645 
0.275 32.841 
0.300 31.622 31.611 31.oll 
0.325 30.107 
0.350 28.470 284470 
0.375 : 26.822 
0.400 252662 25.280 25.269 25.269 
0.425 23.595 
0.450 22.771 22.771 
0.475 21.956 
0.500 21.527 21.504 21.503 
0.525 21.463 
0.550 21.893 21.892 
0.575 22.852 
0.600 26.029 24.477 242426 24.423 
0.625 ‘ 17.574 
0.650 11.927 11.923 
0.675 7.116 
0.700 2.967 2.039 2-832 
0.725 -1.246 
0.750 54404 5.415 
0.775 -9.981 
0.800 0.659 -14.865 -15.245 -154264 
0,825 -21.607 
0.850 -29.337 =29.369 
0.875 -38.893 
0.900 =49.382 -506379 506429 
0.925 63.974 
0.950 -78.942 -79.009 
0.975 -94.170 


1.000 =18.605 106.045 -106.999 107.039 


0.0125 


18.105 
18.727 
20,569 
23,561 
27.594 
32.523 
38.178 
44,375 
50.931 
572668 
64.430 
71.087 
77.538 
63.713 
89.572 
95,100 
100.307 
105.220 
109.8u1 
114.348 
116.688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.488 
146.477 
146.634 
145.972 
144.441 
141.931 
138.265 
133.193 
126.398 
117.500 
106.090 
91.808 
74.474 
54.289 


6.0125 


0.000 
6.197 
i2.let 
17.071 
22-537 
26.606 
29.787 
32.035 
33.371 
33.069 
33.645 
32.441 
31.611 
30.107 
26.470 
26.822 
256269 
23.695 
22.77k 
21.956 
212503 
21.463 
21.892 
22.852 
24.423 
L7.574 
11.923 
7-118 
2-832 
12246 
5.416 
~9.982 
~152266 
-21.609 
-29.371 
-38.896 
-50.433 
-63.978 
-79.013 
-94.174 
-107.042 


0.00625 


16.105 
16.727 
20.569 
23.561 
27.594 
32.523 
386.178 
44,375 
50.931 
57.668 
64.430 
71.087 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109.881 
114.348 
118.688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.488 
146.477 
146.634 
145.972 
144.441 
141.931 
138.264 
133.193 
126.398 
117.499 
106.089 
91.806 
740472 
54.287 


9.00025 


0.000 
6.197 
12.161 
L7.67k 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.841 
3LeGlkl 
30.107 
28.470 
26.822 
25-269 
23.895 
22.77L 
21.956 
21.503 
21.463 
21.892 
22.852 
242423 
17.574 
11.922 
Te1LLd 
2.831 
-1.246 
-5.416 
9.983 
~15.266 
-21.609 
29.372 
-38.8097 
-50.434 
-63.979 
-79.015 
-94.176 
-107.043 


2003125 


16.105 
18.727 
20.569 
23.561 
27.594 
324523 
36.178 
44.375 
50.931 
57.668 
64.430 
TL. O87 
77.538 
83.713 
89.572 
95.100 
100.307 
105.220 
109.4881 
114.348 
118.688 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.487 
146,476 
146.634 
145.971 
144.441 
141.930 
138.263 
133.191 
126.396 
117.497 
106.087 
91.803 
74.469 
54.283 


0003125 


0.000 
6.197 
12.161 
17.671 
22.937 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.641 
3LeOlL 
30.107 
26.470 
26.822 
25.269 
23.095 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.922 
7.118 
2.631 
1.247 
“5.417 
9.983 
-15.267 
=21.611 
-29.373 
38.899 
50.436 
-63.982 
-79.017 
-94.178 
-107.045 


Me 


a a ee bat dat 
ett? ees 


“re 


SECOND ORDER METHOD—~VARIOUS H=--EOQN, (5,46) 
essaBXAMPLE NOs 1) CLEARING TIME 0.60 SECONDS. 


aaaeVALUES OF DELTA ON THIS PADE 


TIME: 


0.00 

0.025 
0.050 
0.075 
06100 
0.125 


206200 


160105 


50.784 


100.316 


136-590 


141.798 


53.011 


0.100 


27.591 


50.920 


17.528 


100.302 


118.687 


136.485 


146.459 


141.892 


117.397 


53.984 


0.050 


4184105 
20.569 
27.594 
38177 
50.930 
64.429 


89.571 


100.307 
409,881 
116.688 
127.304 
136.486 
(63.595 
146.477 
£45.972 
141.932 
133.194 
117.502 

91.808 


54.288 


0.025 


100.307 
105.220 
109.881 
114.348 
119.648 
122.977 
427.304 
131.769 
136.486 
140,662 
143.595 
145.488 
146.477 
146.634 
145.972 
144.442 
141.932 
136.265 
133.194 
126.399 
117.501 
106,091 

91.809 

74.475 

542290 


‘SECOND ORDER METHOC--VARIOUS H=--EQNs. 15.46) 


eee sEXAMPLE NO. 


a eeeVALUES OF OMEGA ON THIS PAGE 


TIME 


0,000 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150. 
0.175 
0.200 
0.225 
0.250 
0.275 
04300 
04325 
0.350 
0.375 
0.400 
0.425 
0450 
0.475 
0.500 
0.525 
0.550 
0.575 
0.600 


0.625 


1.000 


H=0.200 


0.000 


33.504 


25.316 


24.562 


-15.363 


-105.488 


0.100 


0.000 


22.541 


33.376 


31.615 


25.272 


21.507 


244428 


2.822 


-15.319 


-50.53T 


-107.079 


0.050 


0.900 
12.161 
22.537 
29,788 
33.371 
33.645 
31.611 
28.470 
254269 
22.771 
21.503 
21.892 
24-423 
11.923 
2.632 
5.416 
=15.265 
-29.370 
50.429 
-79.007 
-107.039 


1, CLEARING TIME 0.60 SECUNDS. 


0.025 


0.000 
6.197 
12.lol 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.069 
33.645 
32.841 
31.611 
30.107 
28-470 
26.622 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
240423 
17.574 
11.923 
7.118 
2.832 
-1.246 
5.416 
9.982 
-15.265 
-21.608 
-29.370 
-38.895 
504432 
-63.9TT 
-79.012 
-94.173 
-107.041 


0.0125 


18,105 
18.727 
20.569 
23.561 
272594 
32.523 
38,178 
44.375 


89.572 

95.100 
100.307 
105.220 
109,661 
114.348 
118.688 
122.977 
127.304 
131.769 
136.486 
140,662 
143.595 
145.488 
146,477 
146.634 
145.972 
144.441 
141.931 
138.264 
133.193 
126.398 
117.500 
106.090 

91,808 

74.414 

54.289 


0.0125 


0.000 
6.197 
12,161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.841 
31.611 
30.107 
28-470 
26.822 
25.269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
244423 
17.574 
11.922 
7.118 
2.832 
-1.246 
-5-416 
-94982 
15.266 
-21.609 
-29.371 
38.896 
-50.433 
63.978 
=79.013 
-94.174 
- 107.042 


0.00625 


18.105 
186.727 
20,569 
23.561 
276594 
32.523 
36.178 
44,375 
50.931 
57.668 
64.430 
71.087 
77.536 
83.713 
69.572 
95.100 

100.307 

105.220 

109.681 

114.348 

116.688 

122.977 

1272304 

131.769 

136.486 

140,662 

143.595 

145.488 

146.477 

146.634 

145.972 

144.441 

141.931 

138.264 

133.193 

126.398 

117.499 

106.089 
91.806 
74.472 
54.287 


0.00625 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
334645 
32.841 
JLe6Lk 
30.107 
28.470 
26.822 
252269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.922 
7.118 
2.831 
1.246 
52416 
9.983 
~15.266 
-21.609 
-29.372 
-38.697 


50.434 


=-63.979 
-79.015 
94.176 
-107-.043 


+003125: 


95.100 
100.307 
105.220 
109,881 
114.348 
118.668 
122.977 
127.304 
131.769 
136.486 
140.662 
143.595 
145.487 
146.476 
146.634 
145.971 
144.441 
141.930 
138.263 
133.191 
126.396 
117.497 
106.087 

91.803 

74.469 

54.283 


2003125 


0.000 
6.197 
12.161 
17.671 
22.537 
26.608 
29.787 
32.035 
33.371 
33.869 
33.645 
32.841 
31,611 
30.107 
28.470 
26.822 
254269 
23.895 
22.771 
21.956 
21.503 
21.463 
21.892 
22.852 
24.423 
17.574 
11.922 
7.118 
2.831 
-1.247 
-5.417 
-9.983 
-15.267 
-21.611 
-29.373 
-38.899 
-50.436 
-63.982 
-79.017 
-94.178 
-107.045 


NEW AUTOMATIUL PRUCEDURE, 


EXAMPLE 


TIME, 
SEC 


02.000 
0.025 
0.2050 
0.075 
QO.100 
0.125 
0.150 
el TS 
0.200 
04225 
0.2250 
Oe ee 
0.300 
Veaeo 
0-350 
Coats 
024C0 
0.425 
0.450) 
O.475 
0.500 
De De 
0.550 
Ls ous oF 
0.6C9 
0,625 
2650 
O.e075 
OQ. /700 
Dales 
Ge 750 
esta 
0.28C9 
Pe de 
02850 
O0.U7 
0.9CO 
0.925 
0.950 
0.975 
1.000 


NO. 


ly 


VALUES OF 
DELTA 


Th.105 
HH Ht HHH OH 
20.569 
HHH HHH HH 
HAR HG HH Ht 
HE 
38.178 
ere ee 
HH HH aE HE Ht 
Ho Me AE HE Ht 
64.433 
HUH HH RH 
HR HM EH 
He dE He HEA aE at at 
B9.2642 
HRM HE OE 
HH HHH K RF 
HoH He He oe at 
HR RH RH 
KH HH ee 
ce ee 
ae 
Lied g 4 1 
RM HR 
136.461 
Hoe RE 
eee 
HHH HK HH 
146.026 
ee ee ee 
Ho HH GE 
Ho Hea OK Gk 
140.488 
aes. ee 
Howe Ha EE 
ere ee es 
Lis sme 
ae eee ee 
HOM HK Ht HE 
HHH HKG 
46.591 


EQN. 


(6.17) 


CLEARING TIME 0.60 SECONDS. 


VALUES UF 
OMEGA 


0.009 
ees 
Ne ae DEER 
HR HH 
HM HH HH 
ee ee 
29-.(38 
eee ee es 
HOME SEE HE 
Hoot He de Sk Hae 
S3elee 
HOME GE HK 
ROMS He 
Hh Hae ea HE 
206005 
Ke HK Heat 
ee ee es 
re ee 
H St Ht HE HR aE Ot 
ust de He 4h Ht Ht 
HoH te te ot tt 
HR EH RR 
Lilie oS 

Ht te tt He te HH Ht 
23.998 
eH FER Ht 
HK HG HE Se 
oe a 2 
Payer wi 
HR RH H 
$e 4 GE HO HE Ott 
a St te He St oH de 
<=) We RS AY 
HMM te ae eae aE 
a ee eae 
He 4e He He KH 
—-54.84°9 
HHH Ha HE aE 
HH GE Ha HE 
Kod eH de ae at 


wood KE Us 8 by 


AE 


= 


OeLOE 02 


ORAS at at 

QO. 100 
ee ee 
HH KH 
ce 

C.200 
HoH He Ht St at 
He EHH 
HERR HH 
HH HHH 
H HM aE Ht at 
a 
oH tt 

Ge 200 
Hit HOH at 

C.100 
HHH HHH 
Ho ab at ae Ht 
Hea HE HE a Ht 

C.100 
HK HH 
eH eH 
HH HHH 

C.100 
KKH HHH 
HOR Gt aE 
eee 

C.100 
a 
HoH at at at at 
HoH He te 


C.100 


020.0 
oneeee 
oo] «0 
en eaee 
siunee 
saunas 
(0065.0 
eornee 
wenageen 
' pepe 
vos.d 
eehoue 
oseoue 
seeded 
00s: .9 
soanae 


Ce a 


aenene 
axuaee 
nuusee 
aeusraen 
eHarvace 
OUS 66 

sevens 
OG1.9 

anaean 
eeneee 
sonoee 
QOL .0 


Poeones | 


eaeene 
eueene 
OO1.0 

aawasn 
wens e 
PTerie) 
001.9 

ebuten 
weenes 
ebveee 


001.9 


pospuae® 
VERE eee 
ena ce of 


gti e eae 


€euI$ 
enon eee 
BEL < Es 
Hennes eoe 
ronan eee 


NEW AUTOMATIC PROCEDURE, 


EXAMPLE NO. 


TIME, 
SEC 


0.000 
0.025 
0.2050 
0.075 
0.100 
0.125 
0.150 
O-175 
0.2C0 
e225 
G.2250 
Oa is) 
0.300 
soo 
0.350 
OPE fy 
O.4C0 
02425 
0-450 
02475 
0.35C0 
Oe 
0-550 
0.6CO 
e625 
2 650 
0.675 
0.700 
Uel25 
0.750 
ee be pec 
0.8C0 
Oet25 
O.2850 
0.2875 
0.9C0 
02925 
0.2950 
0.975 
1.0C0 


ly 


VALUES CF 


DELTA 


18.105 
HH MRE HEE 

202569 
eH 
HH HH HH 
oe ee 

38.169 
He HR 
eee ee 
HHH HK 
Ke RH KR HK 
He HRM KR EH 
ee 
Ho te 

89.325 
HR HK HT 
sRRURHER 4 
tt HH tt af te dE te 
a ee ee 
HH de HH OH 4 Ht 
HM HR EH 
a 

126.809 
Ho HH HH GE x 

L360. L9%4 
HM OR HK at ae 
Ht tH a HK 
Hod ot de tk RO Ot 

146.447 
ce ee 
he HH K KO 
Hott Mott Hit tt vt 
HOM Ht tt oH Th Ht at 
oe a ee a) 
How RM HAE aE 
He M tt te KO 

LivetSs 
HH HE ME 
eee. ee 
oe ee 


nS eis 3 


EQN. 


(6.218) 


CLEARING TIME 0.60 SECUNDS. 


VALUES OF 


OMEGA 


0.000 
HoH RH GE 
12.161 
HKG Ke He at 
ee ees 
ee a 
29. f9¢ 
oe eee es 
HA a 
ee ee 
HIE ge HEHE aE at 
HE He a aE He 
MoM te hat Kat tt 
oe 
28.346 
ee ee 
RRA K BSE OR 
Mat de A de ae ae Ht 
Heat se a aS a He 
ee ee 
eR aE 
RK Hae eK 
21.602 
ee ee es 
24.2893 
ee ee 
eee 
HK be a 
3.100 
oa aS aS aE at 
Hea a6 Rott a aR 
a eee a 
ce es 
it oH tk ae a a at 
He te ak Ke Hae it 
ee ee 
—“49.564% 
er ae 
OK HE a aE 
Hoa a aE Ht 


~106.231 


A 


E 


0.10E 02 


2 TEP 
S26 


0.050 
ee a 

0.100 
HH HE 
HRS 
HM HE Hat 
GC.200 
HH 
HE EH 
He eR 
a. 
Ree HH 
tH 
ee es 

QO.200 
Kab aE at 
He a aE 
ee ee | 
HK HF 
He 
aed at at 
ee es 

0.200 
tt Me ae ae at 

QO. 100 
ee 
ee ee 
4 ah at at Ht it 

0.200 
Hook at eH 
i it at He at 
Hoe tr oe a 
eae a Ht 
a 
ee 
Me EK 

0.200 
HH tt ae at at 
a 4 th HE at ae 
Hie ie It 


C.200 


020.0 

aeenes 
001.0 

+e eee 
aoxneae 
enaete 
OOS .0 

sou ene 
seunae 
asateoe 


sones 


ensoae 
arene ® 
aieaven 
oe nea 
ree 2 oh 
Pr Se 
+6 eH OM 
onan 
ee 
w#encee 
OHS .O0 

aan eww 
O01. 

auveura 
enecua 
Pe 
ODS .0 

aranTen 
one oe 
auveter 


ao ee 
Pe en 
ota, 
corer ore 
Ce ee 
ssnee ane. 


ceanasey 


DEVOGELAERE METHOD--VARIOUS H-~EQN. Asl-.5 
eeeeEXAMPLE NOs 14 CLEARING TIME 0,60 SECONDS, 


+eeeVALUES OF DELTA ON THIS PAGE 
TIME H=0.200 0.100 0.050 0.025 0.0125 0.00625 .003125 


0.000 18,105 18,105 18.105 18.105 18.105 16.105 18,105 


0.025 16.727 18.727 18.727 18.727 
0.050 20.568 20.569 20.569 20.569 20.569 
0.075 23.561 23.561 23.561 23.961 
0.100 27.569 27.593 27.594 27.594 27.594 27.594 
0.125 32.523 32.523 32.523 32.523 
0.150 38.175 48.178 38.178 38.178 38.178 
0.175 44.375 44.375 44.375 44.375 
0.200 49.463 50.856 50.926 50.930 50.931 50.931 50.931 
0.225 57.667 57.668 57.668 57.668 
0.250 64.425 64.430 64.430 644430 64.430 
0.275 71.067 71.087 71.087 71.087 
0,300 17.442 77.533 77.538 77.538 77.538 77.538 
0.325 83.713 83.713 83.713 83.713 
0.350 89,567 89.571 89.572 69.572 89.572 
0.375 75.100 95.100 95.100 95.100 
0.400 98.010 100.212 100.303 100.307 100.307 100.307 100.307 
0.425 j 405.219 105.220 105.220 105.220 
0.450 109.878 109.881 109.861 109.881 109.881 
0.475 114.348 114.348 114.348 114.348 
0.500 118.579 118.684 118.688 118,688 118.688 118.688 
0.525 (22.977 122.977 122.977 122.977 
0.550 127.300 127.304 127.304 127.304 127.304 
0.575 131.769 131.769 131.769 131.769 
0.600 1324020 136.317 1364480 136.485 136.486 1364486 136.486 
0.625 140.660 140.662 140.662 140.662 
0.650 143.556 143.589 143.594 143.595 143.595 
0.675 145.478 145.486 145.487 145.487 
0.700 145.550 146.363 146.462 146.475 146.476 146.476 
0.725 146.613 146.632 146.634 146.634 
0.750 145.752 145.944 145.968 145.971 145.971 
0.775 144.404 144.437 144.441 144.441 
0.800 107.157 138.606 1416542 141.841 141.925 141.930 141.930 
0.825 138.199 138.256 138.263 138.263 
0.850 132.526 133.108 133.182 133.191 133.191 
0.875 126.288 126.384 126.396 126.396 
0.900 10H.221 116.399 117.359 117.482 117.497 117.497 
0.925 105.913 106.067 106.086 (106.086 
0.950 90.106 91.590 91-740 91.803 91.803 
0.975 7446215 74.441 74.468 74.468 
1.c00 -37.090 36.058 52.010 53.995 54.251 54.282 54.282 


DEVOGELAERE METHOU--VARITUUS H--EQNe. Act-.5 


eeesEXAMPLE AUs Ly CLEARING TIME 0.60 SECONDS. 


eeeeVALUES CF OMEGA UN THIS PAGE 


TIME H=0.2C0 0.1iCU 0.050 0.025 0.0125 0.00625 2003125 
0.c00 0.000 0.000 0.000 0.000 0.900 0.000 0.000 
0.025 Goi 9? Sel di 6.197 6.197 
0.050 12.159 12.160 12086, 12.16 iz.t6l 
0.075 17.67i L7.671 i7.67k 17.671 
0.100 22.498 22.535 22.537 22.537 220537 220537 
0.125 26.608 26.606 2606%; 26.608 
0.150 29.7B5 29.787 29.787 29.787 29.787 
0.175 32.035 32.035 32.035 32.035 
0.200 324334 33.331 33.369 33.371 33.37h 33.371 33.371 
0.225 33.869 33.869 33.609 33.869 
0.250 33.644 33.645 33.645 33.645 334645 
0.275 32.841 32.841 32.841 32.841 
0.300 31.602 31.611 3L.6ll 31.611 31.611 3L.6LL 
0.325 30.107 30.107 30.107 30.107 
0.350 28.470 28.470 28.410 28.470 28.470 
0.375 26.822 26.822 26.822 26.622 
0.400 24,832 25.261 25.269 25.269 25.269 25.269 25.269 
04425 23.895 23.895 23.895 23.895 
0.450 22.771 22.771 22.77k 22.771 22.771 
0.475 214956 21.956 21.956 21.956 
0.500 21.464 21.502 21.503 21.503 21.503 21.503 
0.525 21.463 214463 210463 21.463 
0.550 21.890 21.892 21.892 21.892 21.892 
0.575 22.852 22.852 22.852 22.852 
0.600 21.459 24.314 24.420 24.423 24.423 246423 24.423 
0.625 ’ 17.555 17.572 17.573 17.574 
0.650 11.756 11.903 11.920 11.922 11.922 
0.675 7.096 7ellS 7.118 7.118 
0.700 0.869 2-617 2.805 2.828 2.831 2.831 
0.725 =1.279 -14.250 1.247 1.247 
0.750 -5.746 =-5.458 5.421 5.417 5.417 
BER HB ciesee © -isicer™ Cisszer 
- =19.835 -15.804 -15.335 -15. -15. -15. 
D.825 eb ‘ =21.697 -216620 -21e611 -21.611 
0.850 =30.240 -29.482 -29.385 -294373 29.373 
0.875 ~39.035 -38.914  -38.899 -38.899 
0.900 602979 -510734 -50.600 50.454  -500437 -504436 
0.925 640172 64.003 63.983 634982 
0.950 -80.604 -79.220 -79.040 -79.018 -79.018 


0.975 94.368 -94.199 -94.179 -94.179 
1.000 =68.208 -113.354 108.099 -197.183 -107.060 -107.045 -107.045 


we 6. a 
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a 
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GILL REPEATEO-~VARIOUS H--EQN. (3.23) 


oees EXAMPLE NOw Le CLEARING TIME 0.60 SECONDS. 


TIME, VALUES OF DELTA FOR VARIOUS H 
SEC BEGINNING AT H = 0.200. 

0.000 18.105 18.105 18.105 18.105 
0.025 18.727 
0.050 20.568 20.569 
0.075 23.561 
0.100 27.582 27.592 27.594 
0.125 32.522 
0.150 38.172 38.177 
0.175 44.375 
0.200 50.251 50.819 50.922 50.930 
0.225 57.667 
0.250 64.419 64.429 
0.275 71.086 
0.300 77.348 77.526 77.537 
0.325 83.712 
0.350 39.559 89.571 
0.375 95.099 
0.400 97.598 100.067 100.293 100.306 
0.425 105.219 
0.450 109.866 109.881 
0.475 114.348 
0.500 118.360 118.670 118.687 
0.525 122.976 
0.550 127.283 127.303 
0.575 131.768 
0.600 129.362 135.959 136.459 136.484 
0.625 140.661 
0.650 1432553 143.593 
0.675 145.485 
0.700 145.187 146.416 146.474 
0.725 146.631 
0.750 145.879 145.968 
0.775 144.436 
0.800 100.074 138.820 141.783 141.924 
0.825 138.256 
0.850 132.953 133.182 
0.4875 126.384 
0.900 109.404 117.113 117.482 
0.925 106.068 
0.950 91.215 91.780 
0.975 714.441 
1.000 ~34.415 38.404 53.498 54.252 


eeeeFUR ALLOWABLE ERROR IN CELTA = 1.0 EL.~DEGe.ee. 
THE NECESSARY VALUE OF H WAS 0.050 
THE LARGEST ERRUR ENCOUNTERED WAS 0.755 

THE NUMBER OF FUNCTION ENTRIES WAS 84 
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JOHNSON-WARD METHODC--VARIOUS H--EQN. (5.12) 


eee eEXAMPLE NOw Ly CLEARING TIME 0.60 SECONDS. 


TIME, VALUES OF DELTA FOR VARIOUS H 
SEC BEGINNING AT H = 0.200. 

0.000 18.105 18.105 18.105 18.105 18.105 
0.025 18.727 lag ter 
0.050 20.569 20.569 20.569 
0.075 23.561 23.561 
0.100 27.601 27.595 27.594 27.594 
0.125 32.523 32.523 
0.150 38.182 38.178 38.178 
0.175 44.377 44.376 
0.200 51.315 50.986 50.942 50.932 50.931 
0.225 57.670 57.668 
0.250 64.450 64.433 642430 
0.275 71.090 71.087 
0.300 77.717 77.567 77.542 77.538 
0.325 83.718 83.714 
0.350 89.610 89.577 89.572 
0.375 95.106 95.101 
0.400 102.442 100.62? 100.354 100.313 100.308 
0.425 105.226 105.220 
0.450 109.938 109.889 109.882 
0.47 114.357 114.349 
0.500 119.177 118.759 118.697 118.689 
0.525 122.988 122.978 
0.550 127.395 127.316 127.306 
0.575 131.783 131.771 
0.600 141.646 137.332 136.607 136.502 136.488 
0.625 140.681 140.665 
0.650 143.776 143.619 143.598 
0.675 145.518 145.492 
0.700 148.546 146.781 146.517 1463482 
0.725 146.687 146.641 
0.750 146.498 146.041 145.981 
0.775 144.531 144.453 
0.800 180.694 148.371 142.835 142.049 141.946 
0.825 136.417 136.284 
0.850 134.723 133.391 133.218 
0.875 126.653 126.430 
0.900 136.156 120.012 117.824 117.541 
0.925 106.496 106.141 
0.950 95.682 92.304 91.870 
0.975 75.063 74.548 
1.000 4332669 97.789 59.540 54.956 54.372 


eéeeFOR ALLOWABLE ERROR IN CELTA = 1.0 EL -DEGeceocs 


THE NECESSARY VALUE OF H WAS 0.025 


THE ‘LARGEST ERROR ENCOUNTERED WAS 0.583 
THE NUMBER OF FUNCTION ENTRIES WAS 164 
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CCLLATZ PROCEDURE--VARIOUS H--EQN. (5.17) 


eoeeEXAMPLE NOw Ly CLEARING TIME 0.60 SECONDS. 


TIME, VALUES OF DELTA FOR VARIOUS H 
SEC BEGINNING AT H = 0.200. 
0.000 18.105 18.105 18.105 18.105 
0.025 1as727 
0.050 20.568 20.569 
0.075 23.561 
0.100 27.582 27.593 27.594 
0.125 323523 
0.150 38.175 38.178 
0.175 44.375 
0.200 50.196 50.881 50.927 50.930 
0.225 57.667 
0.250 64.426 64.430 
0.275 71.087 
0.300 77.482 772534 77.538 
0.325 83.713 
0.350 89.568 89.572 
0.375 95.100 
0.400 100.261 100.256 100.304 100.307 
0.425 105.219 
0.450 109.879 109.881 
0.475 114.348 
0.500 118.628 118.685 118.688 
0.525 1222977 
0.550 1275301 127.304 
0.575 131.769 
0.600 136.240 136.386 136.481 136.486 
0.625 140.662 
0.650 143.586 143.595 
0.675 145.487 
0.700 146.143 146.461 146.476 
0.725 146.633 
0.750 145.948 145.971 
0.775 144.440 
0.800 137.499 141.089 141.893 141.930 
0.825 138.263 
0.850 133.131 133.191 
0.875 L2Z6e399 
0.900 115.254 117.399 117.496 
0.925 106.085 
0.950 91.648 91.801 
0.975 74.466 
1.000 30.438 49.436 54.073 54.280 


weeeFOR ALLOWABLE ERROR IN DELTA = 140 EL-DEGeoeee 
THE NECESSARY VALUE OF H WAS 0.050 
THE LARGEST ERROR ENCOUNTERED WAS 0.207 

THE NUMBER OF FUNCTION ENTRIES WAS 63 
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SECOND ORDER METHOC--VARIOUS H--EQN. (5.46) 


eee e EXAMPLE NOw Ly CLEARING TIME 0.60 SECONDS. 


TIME, VALUES OF DELTA FOR VARIOUS H 
SEC BEGINNING AT H = 0.200. 


0.C00 18.105 18.105 
0.025 
0.050 
0.075 
0.100 21,591 
0.125 
0.150 
0.175 
0.290 50.784 50.920 
0.225 
0.250 
0.275 
0.300 77.528 
0.325 
0.350 
0.375 
0.400 100.316 100.302 
0.425 
0.450 
0.475 
0.500 118.687 
0.525 
0.550 
0.575 
0.600 136.590 136.485 
0.625 
0.650 
0.675 
0.700 146.459 
0.725 
0.750 
0.775 
0.800 141.798 141.892 
0.825 
0.850 
0.875 
0.900 117.397 
0.925 
0.950 
0.975 
1.000 53.011 534984 


eeeeFOR ALLOWABLE ERROR IN OELTA = 1.0 EL.DEG..... 
THE NECESSARY VALUE OF H WAS 0.200 
THE LARGEST ERROR ENCOUNTERED WAS 0.973 

THE NUMBER OF FUNCTION ENTRIES WAS 20 
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DEVOGELAERE METHCO--VARIOUS H--EQN. Awl-.5 


eeeeEXAMPLE NO, 


TIME, 
SEG 


0.000 
0.025 
0.050 
0.075 
0.100 
0.125 
0.150 
0.175 
0.200 
0.225 
0.250 
0.275 
0.300 
0.325 
0.350 
0.375 
0-400 
0.425 
0.450 
0.475 
0.525 
0.550 
0.575 
0.600 
0.625 
0.650 
0.675 
0.700 
0.725 
0.750 
0.775 
0.800 
0.825 
0.850 
0.875 
0-900 
0.925 
0.950 
0.975 
1.000 


eeeeFUR ALLOWABLE ERROR 


18.105 


49.463 


98.010 


132.020 


107.157 


ete 3 


CLEARING 


18.105 18.105 
20.568 
21.6569 214593 
38.175 
50.856 50.926 
64.2425 
17.442 77.533 
89.567 
100.212 100.303 
109.878 
118.579 116.684 
127.300 
136.317 136.480 
143.556 
145.550 146.363 
145.752 
138.606 141.542 
132.526 
108.221 116.399 
90-106 
36.058 52-010 
IN DELTA = 


TIME 0.60 SECONDS. 


VALUES OF DELTA FOR VARIOUS H 
BEGINNING AT H = 


0.200. 
18.105 18.105 
18.727 18.727 
20.569 20.569 
23.561 23.561 
27.594 27.594 
IeeJes g£eJ¢3 
38.178 38.178 
44.3175 44.375 
50.930 50.931 
57.667 57.668 
64.430 64.430 
T1.087 T1.087 
77.538 77.538 
83.713 B4e113 
89.571 89.572 
95-100 95.100 

100.307 100.307 
105.219 105.220 
109.881 109.881 
114.348 114.348 
118.688 118.688 
eae ae f 122.977 
127.304 127.304 
131.769 Ag1e769 
136.485 136.486 
140.660 140.662 
143.589 143.594 
145.478 145.486 
146.462 146.475 
146.613 146.632 
145.944 145.968 
144.404 144.437 
141.881 141.925 
138.199 138.256 
133.108 133.182 
126.288 126.384 
ihe OF rye 117.482 
LOS eal 106.067 
91.590 91.780 
74.2215 74.2441 
SEAM E: 54.251 


1.0 EL-DEGseeee 
THE NECESSARY VALUE OF H WAS 0.025 


THE LARGEST ERROR ENCOUNTERED WAS 02256 
THE NUMBER OF FUNCTION ENTRIES WAS 
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